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Abstract 



Four point correlation functions for many electrons at finite temperature in 
periodic lattice of dimension d (> 1) are analyzed by the perturbation theory 
with respect to the coupling constant. The correlation functions are charac- 
terized as a limit of finite dimensional Grassmann integrals. A lower bound on 
the radius of convergence and an upper bound on the perturbation series are 
obtained by evaluating the Taylor expansion of logarithm of the finite dimen- 
sional Grassmann Gaussian integrals. The perturbation series up to second 
order is numerically implemented along with the volume-independent upper 
bounds on the sum of the higher order terms in 2 dimensional case. 
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1 Introduction 

The thermal average of an observable O for many electrons in a solid is expressed 
as Tr e~P H 0/ Tr e~^ H , where H is a Hamiltonian representing the total energy of 
the system, (3 is the inverse temperature and the trace operation Tr is taken over 
the Fermionic Fock space, the Hilbert space of all the possible states of electrons. If 
the movements of electrons are confined in finite lattice sites under periodic bound- 
ary condition, the Fermionic Fock space becomes finite dimensional. The thermal 
average Tr e- (3H 0/ Tr e~ pH is defined as a quotient of finite sums over the orthonor- 
mal basis spanning the space. Though the expectation value Tr e~ /3H Of Tr e~@ H 
has a clear mathematical meaning in this setting, to rigorously control its behavior 
for interacting electrons poses a challenge. The purpose of this paper is to ana- 
lyze the thermal expectation value for 4 point functions modeling paired electrons' 
condensation by means of the perturbation theory. 

In the earlier article |12j Koma and Tasaki rigorously proved upper bounds on 

2 point and 4 point correlation functions for the Hubbard model and concluded the 
decay properties in 1 and 2 dimensional cases. In an abstract general context, on the 
other hand, Feldman, Knorrer and Trubowitz gave a concise representation of the 
Schwinger functionals formulating the correlation functions via Grassmann integral 
and established upper bounds of the Schwinger functionals in [5] . Let us also remark 
the intensive renormalization group study by the same authors in [5] , which analyzes 
the Grassmann integral formulation corresponding to the temperature zero limit of 
the correlation function for the momentum distribution function. The work [8] was 
presented as the 11th paper in the scries of Feldman, Knorrer and Trubowitz's 2-d 
Fermi liquid construction. A flow chart showing the hierarchical relation between 
these 11 papers is found in the digest [7|. 

In this paper we focus on the correlation function Tr e~° H 0/ Tr e~P H for 4 
point functions O = V ; i 1 f'0i 2 j.^y2j.^yiT +'^y 1 f Vy 2 !'0x2|V'xiT and the Hubbard model 
H defined on a finite lattice. We expand the 4 point correlation function as a 
perturbation series with respect to the coupling constant and study the properties 
of the perturbation series. We especially aim at establishing upper bounds on the 
sum of higher order terms of the perturbation series so that one can numerically 
measure the error between the correlation function and the low order terms of the 
perturbation series. More precisely, our goal is set to 
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(1) find a constant r > such that for any UsR with \U\ < r 



Tr(e-P H C 
Tre-^ 



O) 



OO 



where J7 denotes the coupling constant and a n £ M (Vn € MU {0}), and to 

(2) establish an inequality of the form that for any U 6 K with |J7| < r and 
m S NU {0} 



The inequality claimed in @ is proved in Theorem 14.101 as our main result and 
a volume-independent r required in ([TJ)— (J21) is obtained in Proposition 15.11 for 2 
dimensional case. 

Our strategy is based on the discretization of the integrals over the interval of 
temperature appearing in the temperature-ordered perturbation series. By replac- 
ing the integrals by finite Riemann sums we obtain a fully discrete analog of the 
perturbation series in which all the variables run in finite sets. The discretized per- 
turbation series is formulated in a finite dimensional Grassmann Gaussian integral, 
which is rigorously defined as a linear functional on the finite dimensional linear 
space of Grassmann algebras. See [H] for another approach to the finite dimen- 
sional Grassmann integral formulation based on the Lie- Trotter type formula. We 
then rewrite the 4 point correlation function as the Taylor series expansion of log- 
arithm of the Grassmann Gaussian integral. By evaluating the partial derivatives 
of logarithm of the Grassmann Gaussian integral, which were characterized as the 
tree expansion by Salmhofer and Wieczerkowski in |22] , and passing the parameter 
defining the Riemann sum to infinity, we obtain an upper bound on each term of 
the perturbation series of the original correlation function. For completeness of 
the paper and convenience for readers, the derivation of the temperature-ordered 
perturbation series is presented in Appendices. 

As a key lemma we make use of the volume- and temperature- independent upper 
bound on the determinant of the covariance matrix recently established by Pedra 
and Salmhofer in [TB]. Pedra-Salmhofer's determinant bound enables us to find a 
numerical upper bound on the Fermionic perturbation theory in a simple argument. 
As one aim, this paper intends to show a practical application of Pedra-Salmhofer's 
determinant bound. 

Let us note that the lower bound on the radius of convergence of the perturbation 
series proved in Theorem 14.101 and Proposition 15.11 below for 2 dimensional case is 
proportional to (3~ 3 . By applying advanced multi-scale, renormalization techniques 
to the correlation functions of the 2 dimensional Hubbard model, Rivasseau ([19]) 
and Afchain, Magnen and Rivasseau (pQ) proved that a lower bound on the radius 
of convergence is proportional to (log/3)~ 2 , which is larger than our lower bound 
for large f3, i.e, small temperature. In this article, however, we feature calculating 



Tr(e-^ g 0) 
Tre"^ 



m 



$>„E7" <R m+1 (\U\) 



where R m +i(\U\) = 0{\U\ m+1 ) as \U\ \ 0. 



3 



the quantities in a simple manner so that readers can verify the construction of the 
theory by themselves, rather than improving the temperature-dependency of the 
convergence of the perturbation theory via large machinery. 

Our motivation to implement the perturbation theory for many electrons with 
rigorous error estimate numerically was grown amid active research of numerical 
analysis for high temperature superconductivity. The macroscopic behavior of elec- 
tromagnetic fields around a type-II superconductor is governed by a system of non- 
linear Maxwell equations called the macroscopic critical-state models. Prigozhin 
initiated the variational formulation of the Bean critical-state model for type-II 
superconductivity and reported numerical simulations by finite element method in 
|17j . Following Prigozhin's preceding work |17j . finite element approximations of 
various macroscopic models have been studied in rigorous levels up until today. See 
[2], [H] for the latest developments on this subject. In a smaller length scale, the 
density of superconducting charge carriers, the induced magnetic field and motions 
of the quantized vortices in a type-II superconductor under an applied magnetic 
field can be simulated by solving the mesoscopic Ginzburg-Landau models. Numer- 
ical approximation schemes for the Ginzburg-Landau models such as finite element 
method, finite difference method and finite volume method are summarized in the 
review article [3] , which also explains extensions of the Ginzburg-Landau models to 
describe high temperature superconductivity characterized by d-wave pairing sym- 
metry. We now turn our attention to microscopic models governing many electrons 
in a solid and try to approximate the 4 point correlation functions, which are be- 
lieved to exhibit the off-diagonal long-range order as explained by Yang in [24] if 
superconductivity is happening in the system. However, the concept of error esti- 
mate for the numerical computation of the correlation functions formulated in the 
Fermionic Fock space is not yet seen in a mathematical literature as we can see for 
the macroscopic critical-state models and the mesoscopic Ginzburg-Landau models 
today. Hence, in this paper we attempt to propose an error analysis for the nu- 
merical approximation of the correlation functions defined in microscopic quantum 
theory and implement our numerical scheme in practice. 

The contents of this paper are outlined as follows. In Section[5]the model Hamil- 
tonian and the correlation function of our interest are defined. The perturbation 
series of the correlation function is derived. In Section [3] the temperature-ordered 
perturbation series of the partition function is discretized and the discretized par- 
tition function is formulated in a finite dimensional Grassmann Gaussian integral. 
In Section Q] each coefficient of the perturbation series of the correlation function is 
evaluated and upper bounds on the sum over higher order terms are obtained as 
our main result. In Section the perturbation series up to 2nd order is numerically 
implemented together with the error estimates between the 2nd order perturbation 
and the correlation function in 2 dimensional case. In Appendix the standard 
properties of the Fermionic Fock space are reviewed. A self-contained proof for 
the temperature-ordered perturbation series expansion is presented in Appendix IBI 
Finally, the temperature-discrete covariance matrix is diagonalized and its determi- 
nant is calculated in Appendix [Cl 



4 



2 The perturbation theory 



In this section we define the Hamiltonian operator, formulate 4 point correlation 
function governed by the Hamiltonian under finite temperature and expand the 
correlation function as a power series of the coupling constant. To analyze the 
properties of the power series of the 4 point correlation function derived in this 
section is set to be the main purpose of this paper. 

2.1 The Hubbard model 

First of all we define the Hubbard model H as the field Hamiltonian operator on 
the Fermionic Fock space along with various notations and parameters treated in 
this paper. 

The spacial lattice T is defined by T := Z d /(LZ) d , where L(g N) is the length 
of one edge of the rectangular lattice and c?(<G N) stands for the space dimension. 

On any set S we define Kronecker's delta S Xiy (x, y e S) by S x , y := 1 if x 
is identical to y in S, S x , y := otherwise. For example, <5(oo)(ll) = 1 f° r 
(0,0),(L,L)eZ 2 /(iZ) 2 . 

For any proposition A the function 1,4 is defined by 

| 1 if A is true, 

A ' [0 otherwise. 

Using the annihilation operator tp xcr and the creation operator i/i*^, which is the 
adjoint operator of ip xa , at site x g T and spin a e {|, J,}, the free part H and the 
interacting part V of the Hubbard model H are defined as follows. 

Ho:= E n™,yT)r^yr,V:=Uj2r^iMxh (2.1) 

x,yer<T,Te{T,i} xer 

where 

F(xo-,vt) := 5 aiT - tJ2(5 x . y - ej + <5 x , y+ej ) 

V j=l 

d 

— t ■ ld>2 (3x,y-ej-e k + ^x,y-ej+e fc + <5x,y+ej -e fc + <5x,y+ej+e fc ) 

3,k=l 
j<k 

the vectors e T (j e {1, • • • , d}) are given by ej(l) = 5j.i for all j, I £ {1, • • • , d}. 
The parameters t,t',fj,,U <G R are called the nearest neighbor hopping amplitude, 
the next to nearest neighbor hopping amplitude, the chemical potential and the 
coupling constant, respectively. Note that the term representing the next to nearest 
neighbor hopping in F(xa 7 yr) is effective only for d > 2. 
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The Hubbard model H is defined by H :— Hq+V and is a self-adjoint operator on 
the Fermionic Fock space Ff(L 2 (T x {|, |};C)). We summarize the definitions and 
the basic properties of the Fermionic Fock space, the annihilation, creation operators 
in AppendixG] Here we note the fact that dimF/(L 2 (r x {|, J.}; C)) = 2 2Ld < +00, 
which means that any linear operator on Ff(L 2 (T x {|,|};C)) can be considered 
as a matrix. 

Let us prepare some more notations used in this paper. For any linear operator 
A : Ff(L 2 (T x {T,|};Q) -> F f (L 2 (T x {U};C)), Tr A is defined by 

2 2Ld 

TvA:=J2 (4>uMl) Ff > 
1=1 

where (-,-) F is the inner product of F f (L 2 (T x {|,|};C)) (see Appendix EJ and 

{^}?=i is an y orthonormal system of Ff(L 2 (T x {f , J.}; C)). The correlation func- 
tion (A) under the finite temperature T is defined by 

_ Tr(e-P H A) 

where (3 := l/(ksT) > with the Boltzmann constant > 0. 

The momentum lattice T* is defined by T* := (2ttZ / L) d / (2irZ) d . 

For any vectors a, 7 of algebra of length n, let (a, 7) denote Ya=i a (07(0- Let 
(•, -) c „ denote the inner product of C™ defined by (u, v) c „ := Ya=i u (0 v (0 f° r an y 
u,veC". 

For any finite set S, §S stands for the number of elements contained in S. 
Let S n denote the set of all the permutations on n elements for n £ N. 



2.2 The correlation function 

Our goal is to analyze the 4 point correlation function 

( - ^xit' ( Ax2l^ ; y2l' ( /'yiT ^^yi^^/yai^^-l^xiT) by means of the perturbation method with 
respect to the coupling constant U. The correlation function of our interest can be 
derived from the logarithm of the partition function. Let us substitute real parame- 
ters {A Xj y )ZiW } Xj y. ZjWg r(C K) into our Hamiltonian H and define the parameterized 
Hamiltonian H\ by 

H X :=H + V X , V x := ]T f/x, y ,z,w^ T ^^w|^ t , (2.3) 

x,y ,z,w£r 

where we set 

^x,y,z,w • — U S X yS z w^x,z ^x,y.z,w ~t~ ^z,w,x,y; (^-4) 

for all x, y, z,w £ T. Note that H\ still keeps the self-adjoint property and that 

^A|Ax t y,« t w=0,Vx,y,B,wer = H. 

To simplify notations, let X represent a vector in T 4 in our argument unless 
otherwise stated. From now we fix 4 sites Xi , X2 , y 1 , Y2 € T to define the correlation 
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function (^^V^^lVViT + ^l^l^l^l) and write ^ = ( x i: x 2, Yi, Ya) 
and X 2 = (yi,y 2 ,xi,x 2 ). 



Lemma 2.1. The following equality holds. 



1 a 



/3<9A 



log 



•Vi 



Tr e-/ 3 ^' 



\ x =o 



(2.5) 



Remark 2.2. Since H\ is self-adjoint, its spectrum o-(H\) is a subset of M. The 
spectral mapping theorem (see, e.g, Section VIII- 7, Corollary 1]) shows that 
{e~P x } xea (H x ) is the spectrum of e~ f3Hx . Thus, Tr e~ /3Hx > 0. For the same reason 
as above the inequality Tre~P B ° > holds. Therefore, log(Tr e~ rm> - / Tr e _/3H °) is 
well-defined. 

Let C(Ff(L 2 (T x {f, J.}; C))) denote the space of linear operators on Ff(L 2 (T x 
{I, J,}; C)). The proof of Lemma [2~T1 is based on the following lemma. 

Lemma 2.3. Let (a,b) be an interval o/M. Assume that A : (a, b) — > C(Ff(L 2 (Tx 
{t, j};C))) is an operator-valued C l -class function. The following equality holds. 
For all s£ (a, b) 



ds 



A{s) = / e (l-t)A(s)± A{s)e tA( s ) dL 



ds 



Proof. Fix any s £ (a, 6) and take small e > such that [s — e, s + e] C (a, 6). For 
any s' S (s — £, s + e) 



e A(s) _ e A( s ') 



(l-t)A(s) tA(s>)n=l 
e e Jt=o 



1 d 



!L^-t)A{s) e tA{s') )dt 



= / e ^ A ^(A(s) - A(s'))e tA ^dt. 



Moreover, we see that 
d 



ds 



lim 

s'eil-e's + e) 



e A(s) _ e A(s') 



lim 

'€(a-e,a+e) 



1 Jl-t)A(s) A ( s ) ~ A ( s ') tA( s ') 



dt 



s — s' 



where we have used the inequality 
A(s)-A(s') A(sl) 



s — s 







< sup 


as 


6e[s-e,s+e] 



1 e (i-t)A ( s)d tA(s)d ^ 
ds 



with the operator norm | • | and Lebesgue's dominated convergence theorem to 
exchange the order of the limit operation and the integral. □ 
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Proof of Lemma \2.1\ Since the operator- valued function \% i— ► H\ is continuously 
differentiable on any interval containing inside, we can apply Lemma 12.31 to have 



1 d , /Tre"^ 
■log 



(3dX Xi 6 V Tre ~' 3ifo 



A A .=0 



9 



Tr / e^-^-P^^-i-PH^) 



1 \Vo dX Xi 



A A .=0 



/3 Tr e-^ 

i Tr (eCi-tX-/^)^*^*^^ + ^ yiT ^ 2i V'x 2 iV'x 1 T)e t (-' 3H ) 

= (V'^T^xai^i^yiT +V'y lT V'y 2 iV'x 2 iV'x 1 T) J 

where we have used the equality that Tr(AB) = Tr(BA) for any operators A, £?. □ 
2.3 The perturbation series 

The partition function Tr e~@ Hx / Ty e~P Ha can be expanded as a power series of the 
parameter {Ux}x£r 4 - We give the derivation of the temperature-ordered pertur- 
bation series in Appendix IBl Here we only state the result. 

Proposition 2.4. For any [/eR and {Xx}xer 4 C M, 

Tre-fl g * 
Tr e-$ H * 

n=1 ' V X23-l : X23,y2j-l,y2jerCT2j-l,<T2 3 e{T:l} 

• C / x 2j -i,x 2j: y2 3 -i,y2 3 det (C (xj ajXj,ykO-kXk))i<j,k<2n 



x 2j =a: 2j-l 
Vj£{l,2,-. ,n} 



(2.6) 



where the constraint X2j = %2j-i requires the variable xij to take the same value as 
X2j—i for all j € {1, 2, • • • , n} and each component of the covariance matrix 
{C(-x.jO-jXj,y k o-kXk))i<j,k<2n is defined by 

Ci^ax, yry) := ^ £ e ^ e -^ _ ^1-) (2.7) 

with the dispersion relation 

d d 
S k :=-2^cos«k,e i ))-4t'-l d > 2 J] cos«k,e i ))cos((k,e fc ))-/*. (2.8) 

7 = 1 },k=l 
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Lemma 12 . 1 1 indicates that we can construct the power series of 
WiiT^U^l^yiT + ''Pyi&kl^i^*^ b y substituting the series flU} into the 
Taylor series expansion of the function log(a;) around x = 1. Since the radius of 
convergence of the Taylor series of log(x) around 1 is 1, we need to know when the 
inequality | Tr e^ l3Hx / Tr e~P H ° — 1| < 1 holds beforehand. An answer will be given 
to this question in Proposition 12 . 71 below. 

It will be more convenient for our analysis to generalize the problems so that the 
variables U, {\x}xer 4 , {Ux}xer 4 ar e allowed to be complex. We will then recover 
the statements on our original problem by restricting the variables to be real. For 
{Ux}xer 4 C C we dchne a function P{{Ux}x&r 4 ) by the power series of the right 
hand side of (|2.6|) . Let us recall that the real function log(a;) (x > 0) is extended 
to be the complex analytic function log(z) in the domain {z 6 C | \z — 1| < 1} by 
the power series 

°° (I'm— 1 

io g ( z ) = ^L_i_( 2 ,_i)« 

n=l 

In our argument to clarify when the inequality 

\P({Ux }xer*) ~ 1| < 1 (2.9) 

holds as well as in the proofs of other lemmas in this paper, the following lemma 
on the determinant bound on the covariance matrix plays essential roles. 

Lemma 2.5. \18{ Theorem 2.4] For any n € N, (x.j,aj,Xj),(yj,Tj,yj) 6 T x 
U,l}x [0,0) (Vje{l,--- ,n}), 

sup |det((u J -,Vfc> c „C(XjC7-jXj,yfc7> ! j/fc))i<j ) fc<„| < 4", 

Uj , Vj eC" with Hiyllcn, llv^Hcn <1 
Vi€{l, " ,»>} 

where ||u||c« := (u,u)^ 2 for all ueC". 

Remark 2.6. The statement of [HI Theorem 2.4] is on the determinant bound of 
the covariance matrices independent of the spin coordinate. It is, however, straight- 
forward to derive the bound claimed in Lemma [2.5l on our spin-dependent covariance 
matrix from [TH1 Theorem 2.4]. 

We can expand —l//3d/dX^ \og(P({Ux}x£r 4 ))\\ x =o,\/X£r 4 as a power series of 
U as follows. 

Proposition 2.7. Assume that UeC satisfies \U\ < log 2/(16/3L 4d ). Then there 
exists e > such that if {Xx}xGr 4 satisfies \\x\ < £ for all X £ T 4 , the inequality 
holds. Moreover, we have 



i 1 00 

' log(P({U x }xer*)) ^ =0 (2.10) 



(3d\ 
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where the coefficients {a n }$£L are given by 



1 d 



f3d\ Xl 



( n+l 

E 

V i=1 



(-l)i-i 



E 



» fc >l ( Vfee{l,--- ,3} 



n G m 

k=l 



(2.11) 



VArer 4 



with {G n }^ = i defined by 



1 



<>■■ —7 > i E E 

X2j-l:X2j,y2j-l,y2jer<T2j-l : <T23G{T,l} 



n! j=i 



(^X 2 j_l,X 2 j^y2j-l,y2j^X2j-l,y2j-l + A X2 j-_l,X2 3 -,y2j-l,y2j + -^y2j-l,y2j ,X 2 j-l ,X 2 j / 



det(C(xja-ja;j,yfc(T ) fcX)b))i<j i fc<2n 



a3 2j= ;c 2j-l 
Vje{l,2,--- ,»} 



(2.12) 



Proof. Let us fix [/ e C with |C/| < log 2/(16/?L 4d ). Take any 

£ € (0,log2/(32/3L 4d ) - \U\/2) and assume that {\ x }xev* satisfies \\ x \ < e for all 
X e T 4 . Then, we see that for all X e T 4 



1^1 < 



log 2 



16/3L 4d ' 

By using the inequality (|2.13|) and Lemma T2. 51 we observe that 



Joe 2 



(2.13) 



1 = 1. (2.14) 



The inequality (|2.14|) allows us to consider \og(P({Ux}xer 4 )) as an analytic func- 
tion of the multi- variable {Ux}xer i m the domain (|2.13|) . Moreover, we have 



1 d 



pd\ x 



■\og(P({Ux}xer*)) 



Xi 



x x =o 
v;ter 4 



1 °° / 



p 



x x =o 



P({Ux}xer*) 



x x =o 1 
v;ter 4 



(2.15) 



where we used the equality that d\og(z)/dz = X)m=o(~ l) m ( z ~ l) m (V2 € C with 
\z — 1| < 1). Furthermore, we can write 



P({U x }x^) 

d 



x x =o 



1= E G »L=o u n , 

n=l VA-er 4 



(2.16) 



-P({£M* e r 4 ) Xx=0 =E^ G « 

■*i v;ter 4 n= i 



v;ter 4 



IT 
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where G n (n E N) is defined in (|2~T2"ll . By substituting (|^TS|) into (gH]) we obtain 
1 d 



POX* 



■log(P({U x } Xeri )) 



x x =o 
v^er 4 



1 oo / oo 

Ie(-e g » 







IT" 



m=0 



n=l 



v^er 4 



n=l ■*! 



G, 



(2.17) 



A^=0 



71—1 



Again by using Lemma 12.51 we can show that for CeC with 
\U\ < log2/(16/3L 4d ) 



OO 

E 

n=l 









v^er 4 



|C/|"<L 



E 

n=l 



(9 



v^er 4 



< oo . 



(2.18) 



Since the radius of convergence of the power series X)m=o z ' n 1S 1' ^ ne inequalities 
P-18P provide a sufficient condition to reorder the right hand side of (|2.17p (see, 
e -g; [141 Theorem 3.1, Theorem 3.4] for products and compositions of convergent 
power series) to deduce 




n+l—l 

>-x=o E 
vxer 4 j—i 



E 



»iH \-m.j=n+l-l 

n fc >l,V*e{l,- .3} 



Arranging (|2~H?1) yields l[2TTT|) . 



(2.19) 
□ 



By restricting J7 to be real in (|2.10ll . we obtain the power series expansion of 
the correlation function (V'^t^xV'ya-lV'yiT + V' yi T , /'y 2 xV ; x 2 i'0x 1 T}- At this point, 
however, we only know that the series X^^Lo a nU n converges for U G C with \U\ < 
log 2/(16/3L 4d ), which heavily depends on the volume factor L d . With the aim of 
enlarging the radius of convergence and finding upper bounds of the power series 
SJ^Lo a nU n , we will construct our theory in the following sections. 



3 Grassmann Gaussian integral formulation 

In this section we discretize the integrals over [0,(3] contained in the perturbation 
series P{{Ux}xer 4 ) so that the discretized perturbation series can be formulated 
in a Grassmann Gaussian integral involving only finite dimensional Grassmann al- 
gebras. Moreover, by showing that the discrete analog of P uniformly converges to 
the original P, we characterize our partition function Tr er^ H j Tr e~^ H ° and the 4 
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point function (■0xitVx 2 |' ! /V2-IVViT + V , y 1 i-V , y 2 iV , x 3 J.' ! /'xiT) as a limit of finite dimen- 
sional Grassmann integrals. The finite dimensional Grassmann Gaussian integral 
formulation will then enable us to apply the tree formula for the connected part of 
the exponential of Laplacian operator of the Grassmann left derivatives to express 
each term of the discretized perturbation series as a finite sum over trees in Section 

H 

3.1 Discretization of the integral over [0,0] 

We define the fully discrete perturbation series by replacing the integral dx in 
P{{Ux}x£T i ) by the Riemann sum. Let us introduce finite sets [0, (3)h and [— (3, (3)h 
parameterized by h £ N//3 as follows. 

[-p,P)h ■= |-A -/3 + ^ } -^+|,--- ,/3-i 

Note that j}[0, 0) h = (3h and fl[-/3, 0) h = 20h. We define the function P h {{U x }x^) 
of the multi- variable {Ux}xer 4 {<^ C) by 



Ph({Ux}xer±) 

L d f3h 

1 + 

n=1 '"' J ~ \ " X23-l,X23,y2 3 -l,y2iero-2 3 _l,O-2 3 -e{TlJ.}a:2J-l,X2 3 e[0,/9)ft 



L"f3h 1 / 1 

E^n i-- h E E E 



• ^2j-iA S cr 2:j ,iS X2:i _ uX2j U X2j _ uX2jt y 2:i _ u y 2:i j det(C(x.jO-jXj,yk<TkXk))i<j,k<2n- 

(3.1) 

Note that if n > L d f3h, det(C(x.jajXj,ykTkyk))i<j,k<2n = for any (xj,aj,Xj), 
friiWj) e rx{t,I}x[0,^ (j S {1,- •• ,2n}), since ftrx{T,l}x[0,/3)h = 2L d (3h. 

Let us summarize the properties of the function Ph in the same manner as in 
Proposition ^. 71 



Lemma 3.1. Assume that U E C satisfies \U\ < log2/(16/3L 4d ). The following 
statements hold. 

(i) There exists e > such that for any {A;t};tfer 4 (c C) with \Xx\ < £ (VA" S T 4 ) 
and any /i € N//3, £/ie inequality \Ph({Ux} x er 4 ) — 1 1 < 1 holds. 

(ii) For any h £ N//3 

-fl^l°g(iM{tf*}* 6 rO) x x=0 =E a ^ C/n ' 
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where the coefficients {ah,n}^=o are given by 

_ i d i-iy- 1 




(3.2) 



x x =o 



with {Gh, n }^Li defined by 

1 n ( 1 

n\ j=i \ h 



G ^ : =hk\-\. E E E 

X2j-l,X2j,y23-l:y2jer<J 2j -l,(T2jG{T:l} ^2 j - 1 , ?2 j £ [0 , 0) h 



' ^(T2j-l,T^(T2i,i^2;2 3 -l,a;23 

( < 5x2j-l,X2j'5y2j-l,y2j^X2j-l,y2j-l + ^X2j-l,X 2 j,y23-l,y2j "I" -^y2j-l,y 2 j,X2j-l,X 2 j)^ 

• det(C(xLjO-jXj, ykO-kXk))i<j,k<2n- 

(3.3) 

(m,) For all n £ N U {0}, lini; l ^ +00 .; l(E N/ / 3 = a„, where {a„}^L * s defined in 
(ICT1) - (i2~T2ll . 



Proof. The proofs for the claims (JTJ) and {0} are parallel to that of Proposition 12. 71 
based on Lemma [2.51 By the definition (|2.7|) det(C(xj<7j2;j, ykO~kXk))i<j,k<2n is 
piece-wise smooth with respect to the variables {xj}^ 1: which implies that the 
Riemann sums over [0, (3)h in Gh, n all converge to the corresponding integrals in 
G„ as h — ► +oo. Thus, the claim (|Ju|) is true. □ 

Lemma l3~l"l(|m|) tells us that establishing an /i-dependent upper bound on \ah, n \ 
and showing that the upper bound converges as h — » +oo lead to finding a bound 
on \a n \. This goal will be achieved in Section [H 

The main aim of this section is to formulate Ph as a finite dimensional Grassmann 
Gaussian integral, which will be used in the characterization of the coefficients 
Wh^J^Lo in Section[U Though it is not directly required in our search for the upper 
bound on X)^°=o a nU n , to represent the original partition function P and the 4 point 
function (V'xiT' ! /x2.l' i /'y24' ! /ViT +Vy 1 tVy 2 x' i /'x2lV ; xiT) as a limit of the finite dimensional 
Grassmann integrals also interests us. The following uniform convergence property 
of Ph provides a framework to this purpose. The following proposition will be 
referred in the proof of our main theorem Theorem 14.101 as well. 

Proposition 3.2. For any r > 

lim sup \Ph({Ux}xer*)-P({U x } x& *)\ = 0. (3.4) 



U^€C with \U x \<r 



Remark 3.3. For the same reason as for the convergence property Lemma 13.11 
each term of the series Ph({Ux}xer 4 ) converges to the corresponding term of 
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P({Ux}xer 4 ) as h — > +00. By using this fact and Lebesgue's dominated conver- 
gence theorem for ^-space, the convergence property (|3.4|) can be shown. Below 
we present an elementary proof without employing the convergence theorem of the 
Lebesgue integration theory. 

Proof of Proposition By using Lemma 12.51 and the inequality that | Ux | < 
r {VX G T 4 ), we have 



\P({U x } Xer *)-Pk({U x } Xer .)\ < £ ( r/3 L 4d )"4 2n 

* — * rt 1 



n=f3h+l 



El n 
~i n 
„ n\ j=i 



^C 2 j-l,T ( 5cr2j,ll^X2 ; ,-l,X2j,y2j-l,y2j I 
1 X 2 j - 1 ,X 2 j ,y 2 3 - 1 ,y 2 3 G r cr 2 j _ 1 , 0-2 j e { T , 1 } 



s 2j= = 2j-l 

Vie{i,— ,»} 



4 U S det(C( Xj a 

J \ ^2 J -ie[o,/3)h/ 



jXj,yk°'k x k))l<j,k<2n 



a: 2j= x 2j-l 
Vj<={l,-- ,n} 



n=f3h+\ 



0h ~ l r(hj-i + l)/h 



Y,-^r sup fi E / 



ds 2 



3-1 



det(C(x ; ,CT i s : ,-,yfccrfcS fc ))i< :)ifc <2„ 
- det(C(x.jO-jlj/h,y k o-kh/h))i<j,k<2n 



s 2 3 '=»2j-l 

vje{i,- - ,71} 



(3.5) 



We especially need to show that the second term of the right hand side of the 
inequality (|3.5|) converges to as h — ► +00. 

Let us fix n G {2,3,- •• ,(3h} and x J; y 3 G L, ^ G {t,4} (Vj G {1, • • • ,2n}). 
There exists a function 

5 : (-A^"- 1 )/ 2 - M, 5 G C oo (((-/3,/3)\{0})"("- 1 )/ 2 ) 

such that for all s 2j -i G [0, /?) (Vj G {1, • • • , n}) 

ff(si - S3, Si — S 5 , ■ ■ ■ , Si - S 2n -1, S3 - S 5 , • • • , S3 - S 2 „-l, • • • , S 2 „_ 3 - s 2n -i) 
= det (C(XjCTjSj, yfc<TfcSfc))l<j,fc<2n| s 2 3 = = 2 3 -i . 

Note that by using the property that E^ = E^ for all k G T*, we can show 
C(xax,yTy) € R for all (x, a, x), (y,T,y) G L x {f, |} x [0, /3)ft. Thus, the function 
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g is chosen to be real- valued. Then we see that 

■ | det(C(-x.jajSj,yk<TkSk))i<j,k<2n\ V S . 2 | { = 1 S2;, '- 1 } 
- det(C(x i cr i l i //i,y fe cr fc ; fc //i))i< i 

Vje{l, -- ,n} 

22 / ds2 ^i (xi + (i - xi)xi, s ) 

2j _ 1 =o"'^-i/' 1 y 

• - s 3 , • • • ,s 2 „_3 - s 2n -i) - g(h/h- h/h, ■ ■ ■ ,l 2 n-3/h- hn-l/h)\, 

(3.6) 

where the functions xi, Xi,s arc defined by 




Xi := 
and 



J 1 if there exist j, k £ {1, • • • ,n} such that j 7^ fc and ?2j-i = ^2fc-i, 
otherwise, 



1 if s 2 j-i - s 2 fe_i 7^ hj-i/h - hk-i/h 
Xl,s := <( for all j, fc G {1, • • • , n} with j ^ k, 

otherwise. 

Let us fix 1 = (h,h, • • • , Z 2 „-i) and s = (s i; s 3 , • • • , s 2 „-i) with 
fcj-i € {0,1,--- 1}, s 2 j_i € {hj-i/h,(l 2j -i + l)/h) for all j e {l,--- ,n} 

satisfying xi = and Xi.s = 1- In this case l2j-i 7^ ^2fc-i and S23-1 — S2k-i 
hj-i/h — hk-i/h for ah j,k G {l,--- ,«} with j 7^ fc. Note that if hj-i < 
hk-i, hj-i/h — hk-i/h, S2j-\ — S2fe-i € (— /?, 0). If hj-i > hk-i, hj-i/h — 
hk-i/h, s 2 j-i — s 2 /t-i € (0, /?). Let us set the interval I{b, c) for 6, c e M with 7^ c 

by 

7(6, c) := [6, c] if b < c, [c, 6] if b > c. 

Then we see that J(s2j-i — S2k-i,hj-i/h — hk-i/h) C (— (3, /?)\{0} for all j, k £ 
{l,-- - ,n} with j ^ k. Since .9 e C°° (((-/?, /3)\{0})™ (rl_1)/2 ), the mean value 
theorem ensures that for any j, k € {1, • • • , n} with j < k there exists ^ 2 j-i,2fe-i € 
H s 2j-i - S2k-i,hj-i/h — hk-i/h) such that 

g(si -s 3 ,--- ,s 2n -3 -s 2 „-i) - g(h/h- h/h,- ■ ■ , hn-a/h - hn-l/h) 

= (V<7(01,3, • • • , 02n-3,2n-l), 

(si - s 3 - {h/h - h/h), • • • , s 2 „- 3 - s 2 „-i - {hn-a/h- hn-i/h)) 1 ), 
which leads to 

\g(si - s 3 , • • • , s 2 „-3 - s 2n -i) - g(h/h -h/h,-- - , hn-a/h - hn-\/h)\ 

1 /n(n-l)\ 1/2 .„ , „ (3.7) 

<t 9 - sup |V 5 (s)|. 
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Moreover, by using Lemma 1231 we see that for j < k, 

1 1 







9(s2j-l - S 2 fe-l) 

d 



^ E E 

Pi — 0p2— 



9(s 2j -i - s 2 fc_i) 
5 



d(s 2 j-i - s 2 k-i 



g(s\ - S 3 , • • • , S 2n _3 - «2n-l 

+pif2j-l+ Pl S2j-l, y2fe-l+p 2 fJ 2fc-l+p 2 S 2 fc-l) 
+pi f 2fe — 1+pi S2k— 1 j Y2j — 1+P2 1+P2 ^2j — 1 



< 8 • 4 2 ™- 1 sup 



xer,a;e(-/3,/3)\{0} 



<9l 



C(xTx,OT0) 



l2n-l 



(3.8) 



By (|37f [l - (;3T8jl we have 

/ .(i 2J _ 1 +l)/ft 



rfs 2j -i (1 - xi)Xi,< 



n E 

• |5( s i - S3, • • • j s 2n -3 ~ s 2n ~i) - g(h/h- h/h, • • • , hns/h - l 2n -i/h)\ 



< -i sup 

« xer,xe(-/3,/3)\{o} 

On the other hand, note that 



^C(xtx,OT0) 



(Bh) n - 



Ui/h e [0,(3)1 I xi = i} = W - tt{l/^ e [o,/3)£ I xi = 0} 

0h 

n 



where we used the inequality 



(3.10) 



N 
n 



i\ < n 2 N n -\ 



which holds for all N e N and n e {0, 1, • • • , N}. By using Lemma 1531 and (|3~TU|) 
we obtain 



n E 



_ dS2j-l | XI 

\g(si - s 3 ,--- ,s 2n _ 3 - s 2n _i) - g(h/h- l 3 /h, ■ ■ ■ ,l 2n - 3 /h- l 2n -i/h)\ 



< 2h~ n n 2 {l3h) n - 1 4 2n = -n 2 f3 n -H 2n . 



(3.11) 
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Combining ([3j|, $B%, (|3"TTj) with (33]) shows 

sup \P({Ux}xa+)-Ph({Ux}xG«)\ 

U X £C with \U X \<T 

v^er 4 

oo Ph 

< y 4(# u )"^" + -r^('-i 4d r 

z — ' ft! ft ' n! 

n=2 



n=/37i+l 

• ( n(n- l)/3™4 2 ™ sup 
\ xer,xe(-/3, / 3)\{o} 

->o, 

as ft — > +00, ft € N//3. 
Corollary 3.4. for allU eR 



^C(xtx,OT0) 



• 2n a /3 n-1 4 an 



□ 



1 g^-PhC^A^Afer*) 
- — lim 



A^=0 

(3.12) 



Proof. The relation (12. 4[) and Cauchy's integral formula ensure that for any 
{Ux}xev i C C and r > 



0A 



(P^ - P)({£M* 



a 5 



er 4 J 



+ 7^77" (Ph-P){{Ux}x^) 



Xx 
1 

2^i 



dU 



U Xl -U Xl \=r 



A', 



(P^-PH^W^ 



(3.13) 



v^er 4 



By applying Proposition 13.21 to (|3.13p we can show that for any f > and any 
{Uxjxer* with \U X \ < f (V-Y e T 4 ) 



(9 



c9 A 



(P^-P)^}^) 



< - SUP |Pfc({EM*er«) ~ P({Ux}x^)\ -> 

as ft — > +00, ft e N//9. Combining (I3.14|) with Lemma |2~T1 yields (I3.12|) . 



(3.14) 



□ 
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3.2 The Grassmann Gaussian integral 

To deal with the discretized partition function Ph rather than P is advantageous 
since the variables run in the finite set T X {}, }} X [0, (3)h in every term of the power- 
series Ph- Accordingly, we can formulate P% as a Grassmann Gaussian integral on 
finite Grassmann algebras. Elementary calculus on finite Grassmann algebras has 
been summarized in the books [6] , [20] . For a convenience of calculation, especially 
in order to refer to Proposition ! C . 71 shown in Appendix [Cl we assume that h G 2N/ /3 
from now. 

Let us number elements of the set T x {1,4} x [0,/3)/i so that we can write 
rx {4,4} x [Q,(]) h = {(x^a^Xj) | j e {1,- •• ,N}} with N := 2L d f3h. We then in- 
troduce a set of Grassmann algebras denoted by {ip Xj a jXj , ^ x -<t x- \j e {1, • • • , N}}. 
Remind us that the Grassmann algebra {^)yc jajXj , "4> x a -x I J "= {!> ' ' ' > N}} satisfies 
the anti-commutation relations 

for all j,k € {1,---_,JV}. 

Let Cfi/'xjo-jXj , ^x-ff-a- b e {1j * " j A 7 }] denote the complex linear space spanned 
by all the monomials consisting of {ipK f <T ) x :j ,'>P X a-x \j e {!>'" ,^}}- As a linear 
functional on <C[ip XjajXj , ip XjCTjXj \ j € {1, • • • , N}], the Grassmann integral 
/ - d ^ XN a N x N ■ ■ ■ #x l0 ^ 1 #x Jv a J vx JV ■ ■ ■ dipwiX! is defined as follows. 

J V'xiffiXi ' ' ' 4'x N a N x N ' l P Xl cr 1 x 1 ' ' ' ^xnctnXn 

■ d ^x N a N x N ■ ■ ■ d ^ Xl<TlXl dljj XN a N x N ■ ■ ■ #x lCT1 X! 
J ^ Xjl a H x H ■ ■ ■ ^*j n <y ]n x ]n ^Px kl cr kl x kl ' ' '^x km a km x km 

' di/> XN a- NXN • • ■ dip Xl(TlXl d'ip XNCrNXN ■ ■ ■ difj Xl(7lXl 

if n 7^ N or in ^ AT, and linearly extended onto the whole space. 
Let us simply write the vectors of the Grassmann algebras 

= (V'xio-ixi, • • • , Tpx N cr N x N ), i>X — vPxktixiI ■ ' ' i *Px N <J N x N )- 

In order to indicate the dependency on the parameter h, we write the covariance 
matrix as 

Ch ■= {C{^-ja X.j,^.k^kXk))l<j,k<N 

and define a 2N x 2N skew symmetric matrix Ch by 




' l Px j a- j x j ' l Px k a- k x k — ' l Px k a- k x k ' l Px j a-jX j , ' l Px j a- j x j ' l P Xk a kXk 
V'xjO-ja;,- V'xfcCTfcXfc ~ ~~ ^ Xk G k x k rxj (T,-ij 
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The diagonalization of Ch is presented in Appendix [C] Here we note the fact that 
det Ch 7^ proved in Proposition IC . 71 to see that Ch is invertible. 

For any ftyxrfx) e t^W,,^ |j £ {1, • • • ,iV}], e /«*Ac) is defined 

by 



e /(^,^) :=e /(o, 0) (^l ( /(^ X! ^)-/(o,o)r) , 



\n=0 



where /(0,0) denotes the constant part of /(ipXi^x)- 
Let us also write in short 

dip x = d^ N a N x N ■ ■ ■ dip xl<TlXl , dtpx = dip XNaNXN ■ ■ ■ dtp Xl(TlXl . 

Definition 3.5. As a linear functional on < C[ip Xj(TjXj ,ip Xja . jXj \j € {1, • • • , N}], the 
Grassmann Gaussian integral J -d^c'h (ipx, ?Px) i s defined by 

(3.15) 

for all f(if>x,il> x ) E C[^ Xj , ^ ajxj |j 6 {1, • • • , N}}. 

Remark 3.6. The denominator of (|3.15|) is non-zero. In fact, a direct calculation 
and the assumption h € 2N//3 show 

-i<(^ A .,^) t ,c- 1 (v.x,? x ) t )^^ x = (detCfc)-^-!)^- 1 )/ 2 = (detCV)" 1 , 



which takes a positive value independent of h by Proposition IC.7I 

The Grassmann Gaussian integral representation of Ph is as follows. 
Proposition 3.7. Assume that {Ux}xeT i (<^- C) satisfies the equality that for all 

x,y, z,w e r 

^x,y,z,w — ^z,w,x,y- (3.16) 

T/ie following equality holds. 



where 

Proof. By substituting the equalities J \d^,c h (ipx, tpx ) = 1 an d 

= det (Cft, (x Jp CT ip X ip , X fcg (7k, )) l<p,g<n 



19 



for any ji, ■ ■ ■ ,j n , ki, ■ ■ ■ ,k n £ {1, 2, • • • , N} (see Problem 1.13]) into (|3.ip . we 
have 

p,({^W) = i+E^,n \-- h E E 

^ ] <5o-2 3 --l,T^cr2j,l^a;2j-l,a;23^X23-l,X23,y2j-l,y2j J 

Z2j-i,Z2je[0,/3)h / 

V'X2„ ( T2„^2„ • • • V , x 1 <ri! Bl ^y 1<7ia . l ' ' ' *Py 2n a 2n x 2n d HC h (i>x3x) 

( Ldph 1 / 1 

h+E^nf-^ E E E 

V 11=1 \ X2 J _l,Xa^,y2i-l,y2fer<7 2 f-l,ff2je{T.4-}i<!2i-l,Xai6[0,/9) l , 

' fi<T2j-l A^<l 2 j ,l&X 2 j-l, X 2 j C^X 2 j-l,X 2 j ,y2j-l,y2j 

• ^Px 2j a 2j x 2j V'xaj-KTaj -1X2^-1 ^yjy-K^- 1X23-1 ^y 2 jcr 2j x 2 j^ J i^Xi^x) 

(3.17) 

To obtain the last equality of (|3.17p we used the equality (|3 . 1 6|) . □ 

As a corollary, our original partition functions and the correlation function are 
represented as a limit of the finite dimensional Grassmann integrals. 

Corollary 3.8. For any (7£l and {\x}xer L C K, the following equalities hold. 

rp —PH\ r _ 

™ pHo = / e^.-.-er^. J r,..-V hlWl w(*X,^) d/1Ch(V , A . ) ^ jc)i (3 . 18) 

h<=2N/(3 

= (JjfL , flfc E I (^^x^ X2 lx^y 2 lx^ yi U +^ yi ^x^y 2 lx^ 2 lx^x) 

(3.20) 
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Proof. Since the relation (|2.4p implies the condition (|3.16|) . we can apply Proposition 
321 and Proposition O to deduce IpTIS]) . The equality (I3J9]) is (f3~T8l) for A* = 
(VAf e T 4 ). Note the fact that 

9 



\ x =0 



= J! (V'x 1 T : i ; V , x 2 _lxV'y 2 i^^yiTx + V'y 1 TxV ) y 2 ^V'x 2 ixV'x 1 Tx) ( 3 - 21 ) 

xe[o,0) h 

where the differential operator d/dX^ i is defined to act on every coefficient of Grass- 
mann monomials in the expansion e^^^ iUxVh - x ^ x '^ xS> (see Problem 1.3]). 
Moreover, by expanding e^xer^ u x v htX (i> x ,ip x ) one can verify the equality 

8 - (3 ' 22) 



The equality (|3~2TJ|) follows from CorollaryEH PropositionOand (|3~2l> (|3~22]) . □ 



4 Upper bound on the perturbation series 

In this section we calculate upper bounds on our perturbation series X^^Lo a nU n by 
evaluating the tree formula for the connected part of the exponential of Laplacian 
operator. In order to employ the Grassmann Gaussian integral formulation of Ph 
developed in Section [321 we assume that h G 2N//3 throughout this section. 



4.1 The connected part of the exponential of Laplacian op- 
erator 

Our approach to find an upper bound on \a n \ of our perturbation series^^ a n U n 
is based on the characterization of the connected part of the exponential of the 
Laplacian operator of Grassmann left derivatives reported in [22j . Let us construct 
our argument step by step to reveal the structure of the problem. 

The Grassmann integral J e^*e r4 UxVh ' x ^ x ^ x ^ d^Ch^x, ipx) can be viewed 
as an analytic function of the multi- variable {Ux}xer 4 - Since 



= 1, 



u x =o 



if \Ux\ is sufficiently small for all X £ T 4 , the inequality 



< 1 
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holds. Thus, we can define a function Wh({Ux} xer 4 ) by 

W h ({U x } xer *) := log ( [ e^^r* u * v ^^^dvcMx^x) 



j id 

which is analytic in a neighborhood of in C . 

Lemma 4.1. For all h G 2N//3 and weNU {0}, i/ie following equality holds. 



0>h,n 



- y 



d n+1 W h 



fi n - i \ t/l> xi,X2,yi,y2 t -"^zi,zi,zi,zi 



a/7 • • • a/7 

t^L'Zi .Zi .Zi .Zi Z n ,Z n ,Z„ ,Z T1 

d n+1 W h 



'®Uyi ,y2 ,xi ,X2 ,zi,zi,zi ' ' ' <9f^i 



Zf. ,z« ,z ? , ,z„ 



[7^=0 



[7^=0 

v^er 4 ' 



(4.1) 



where au,n was defined in (|3.2p - (j3.3p . 

Proof. The Taylor expansion of M / / l ({J7^}^ gr 4) around is given by 



* — * n 



n< 

n=l Xi 



E A dU Xl ■ ■ ■ dU Xn 



0-^=0 
v;ter 4 



U Xl ---U Xn . 
(4.2) 

Fix any [7 G C with \U\ < log2/(16/JL 4d ). Let e > be the constant claimed 
in Lemma [34] (0) . By using a parameter {Xxjxer 4 ^ C) with |A^| < £ (VAf e 
r 4 ) we define the variable {Ux}xer 4 by the equality (12. 4|) . Then the inequality 
|-f\({£ yr ;r};f6r 4 ) — 1| < 1 holds by Lemma l3~Tl fif and the condition (|3.16p is satisfied. 
Thus, Proposition 13.71 ensures that Wh{{Ux} xev 4 ) = \°g(Ph{{Ux}xer 4 )) I0r this 
{Ux}xer 4 - Moreover, by the equalities ()2.4p and (|4.2p we have that 



1 a 



log(Pfe) 



i a 



A A .=0 

v^er 4 



a 



v;ter 4 



1 ^ [/"-I 

7^ (n-1)! 

n— 1 x zi , 



E 



ai^ 



,z„_ier 



X\ ^^zi ,zi,zi.zi ' ' ' $U Zn l : z n _ 1 ,z 7l _i,z n _i 



1/^=0 



9Up^ 2 dU Zl ,zi ,zi ,zi 



• a/7 



z,i-i,z„_i,z„_i,z„_i 



By uniqueness of Taylor series and Lemma 1341 (JTTJ) we obtain (|4.ip . 



□ 



A message from Lemma HU is that upper bounds on \ah, n \ can be obtained by 
characterizing the partial derivatives of Wh at Ux — (\/X e T 4 ), which is the 
way we follow from now. Since \ah,o\ can be evaluated directly from (|3.2p by using 
Lemma [2.51 let us study the equality (|4.ip for n > 1. Fix any n > 1 and use the 
simplified notations defined as follows. 



Z j := ( Zj , z j} Xj , Zj) G T 4 for Zj G T (Vj G {1, • • • , n}), Z := £ G T 4 



(4.3) 
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Set N n+ i := {0, 1, • • • , n}. By noting that 
d 



jeQ dU z . 



n 



e Zxer*U x V h , x (^,iP x ) dfiCh{i!x ^ x) 



u x =o 



d 



e<3 dU Zj J i6N„+ 
for any Q C N n+ i, we see that 

a 



n {l + U Zj V h} zMx^ x ))dix Ch {^x,^x) 



U Z . =0 



jell+i dU z 



n 



■w h 
d 



u x = o 



.. -log / n {i + u Zj v h<Zj (i> x ,i>x)Wc h (i>x,^x) 



—0 



n 



9 



ieN„ +1 ac/ 2j 
/ 

• log 



i+ V / n v h , Zj (ip x ,i>xWc h {iPx,iPx) n 



u z . =o 



(4.4) 



The Grassmann Gaussian integral contained in the right hand side of (|4.4[) can be 
rewritten as follows. 

Lemma 4.2. Introduce Grassmann algebras {'4'x j <T j x j T' l P'L j <T j x j \j G {1; ' ' ' >-^}} 
dexed by q £ N„+i and write 



), i>X = (^xi (7\x\ , ■ ■ ■ ) V'xjv^Tjv ) 



"7"9 



/or a/Z q G N n +i. Let d/dipx> 9/dip x oe ^ e vectors of left derivatives associated 
with 4>x' ^X' respectively. Then, the following equality holds. For all Q C N„+i 



q&Q q£Q 



(4.5) 



where the Laplacian operator A and its exponential e A are defined by 



A :-- 



E 



p,geN„ 



04) ' C M^ 



A' 



2AT(n+l) 

■ a - y -a' 



=0 
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Remark 4.3. When we introduce another set of Grassmann algebras, let us think 
that the complex linear space spanned by monomials of all the Grassmann algebras 
introduced up to this point is defined on the assumption of multiplication satis- 
fying the anti-commutation relations between these algebras. The notion of the 
Grassmann integral J -dipxdipx is naturally extended to be a linear map from the 
enlarged linear space of all the algebras to the subspace without ipx , ipx ■ 
For a monomial <f>j 1 ■ 

(d/d^fa—fa (Ze{l 
d , 



bj n of Grassmann algebras the left derivative 

• , m}) is defined by 







if there uniquely exists 
k £ {!,■■■ ,n} s.t. / =j k , 
otherwise. 



Then, the left derivative d/d<pi is extended to be a linear map on the linear space 
of monomials of the algebras {4>l\iL\- 

The concepts of Grassmann integrals and left derivatives are generally defined 
as operators on Grassmann algebra with coefficients in a superalgebra (see [SJ 
Chapter I]). 

Proof of Lemma \4-.2\ We define another Grassmann algebra {tIx^x} indexed by 
the sets T x {|, J.} x [0,P)h and N„ + i and the associated left derivative 
{d/di] x ,d/drf x } in the same way as {ipx'^x} an d {9/dtp x ,d/dtp x }. Then, we 
see that for any subset Q C N„ + i with Q ^ 



IT Vh,z q {ipx,ip x ) d ^c h (ipx,tlJx) 



n v KZc 



n V KZq 



n V KZq 

96Q 



f_d d_ 

ydrjx' dff x 

d d 



dr fx ' dr lx 
_d d_ 



n X =V X =° 

VqGQ 



n x =V x =o 

VqGQ 



(4.6) 



n\=n- x =o 

V<jgQ 



qeQ 



Vq£Q 



where we have used the equality that 



(see [SJ Problem 1.13]). To verify the equalities (|4.6|) in more detail, see the books 
[E], [201 for the properties of left derivatives. □ 
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By combining (14. 4|) with (|4.5|) wc obtain 
5 



£7^=0 



ieN„ +1 dU z . 



■log 



1+ E e A n/ M ,(4,^) ^.A^ 



V 



5c»„ +1 



C7 Z =0 
^q 

(4.7) 



In order to characterize the right hand side of (|4.7p , let us review the general theory 
developed in [22] by translating in our setting. Consider a map a from the power 
set P(N„+i) of N n+ i to C defined by a(0) := 1 and for Q e P(N„+i)\{0} 



a(Q) :=e A n V h , Zq ^ q x ^ q x ) 
<?eQ 



i> q x =-4> q x =n ' 



By [22l Lemma 1] there uniquely exists a map a c : P(N n +i) — * C such that for all 
QeP(N n+1 )\{0} _^ 

a(Q)= a c(QoMQ\Qo), 



Q cQ 



where minQ stands for the smallest number contained in Q. In [22J, Lemma 2] it 
was proved that the right hand side of (|4.7[) is equal to a c (N„ + i), which is called 
the connected part of the operator e A . The formula for a c (N n+ i) was given in [2"2l 
Theorem 31. We summarize the result below. 



Lemma 4.4. 

n 



Theorem 3] The following equality holds. 



-w h 



u x =o 



ieN„ +1 dU Zj 

]T B (A q , q ,+A q ,, q ) [ ds Yl m*,s)e A{M{T ^ s)) 



TGT(N„+i) 



n K,^(V4^x) 



7reS„+i(T) 



(4.8) 



where T(N„+i) is i/ie se< o/ all the trees (connected graphs without loop) on N„+i, 

t 



(9 



,0. 
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s := (si,-- - ,s n ), S n +i{T) is a subset of S n +i depending on T, <p(T,ir,s) is a 
real-valued non-negative function of s depending on T and ir with the property that 

f ds V #T,7r,s) = l, (4.9) 

.GS„ +1 (T) 

M(T, 7r,s) is an (n+ 1) x (n + 1) real symmetric non-negative matrix depending on 
T , tt, s satisfying M(T, tt, s) 9i<3 = 1 for all q G N„+i and the operator A(M(T, tt, s)) 
is defined by 

A(M(T ;7 r,s)) := £ M(T, tt, s) M A M . 

In order to bound |a^ )Tl |, the tree formula (|4.8|) will be evaluated in the rest of 
this section. 



4.2 Evaluation of upper bounds 

Here we evaluate the tree expansion given in Lemma 14.41 Let us first prepare 
some necessary tools. The following lemma essentially uses the determinant bound 
Lemma 12.51 



< 4 1 



Lemma 4.5. For any I G N and any p m , q m G N n +i, (x Jm , <Xj m , Xj m ), 
(x km ,o- km ,x km ) GTx {T4} x [O,0) h (VtoG {!,••• ,1}) 



,A(Af(T,7r,s))^ 1 ^ 



Proof. Since M(T, tt, s) is a non-negative real symmetric matrix, there are constants 
lq > (g G N„ + i) and projection matrices P g (g G N n +i) satisfying that Pp-Pg = 
for all p,q E N„+i with p ^ q such that M(T, tt, s) = X)g=o Tq-^V Define an 
(n + 1) x (n + 1) real matrix M by M := J2q=o y/lqPq- Then we see that 



M(T,ir,s) = M l M. 



(4.10) 



By writing M = (vo, • • • ,v„) with vectors v q G R n+1 (q G N„+i), the equality 
(I4.10P implies that M(T, tt, s) Pjg = (v p ,v g ) for all p,q G N n +i. The property that 
M(T, vr,s) 9:9 = 1 (Vg G N n+i ) ensures that for all g G N„+i 

|v,| = 1. (4.11) 



Then we observe that 



A(M(T, 1 r,s)), /) p 



~r9i 



#>i ... ^Pi J," 



> 91 



= (-l)^- 1 )/ 2 det((v Ps ,v gt )Cft(x ;;s cr is a;j s ,Xfe t crfe t Xfc t ))i< S)t < ; . 



CT fcm X k„ 



(4.12) 



2G 



By noting (|4.1ip we can apply Lemma l2~5l to ()4.12[) to deduce the desired inequality. 

□ 

One point we need to carefully deal with in the evaluation of the right hand side 
of (|4.8p is the combinatorial factor, which comes in the expansion of 
H{q,q>}eT(A q , q > + A q i i9 ) II 9 eN„ + i Vh,z q (ipx,4> q x)- ln order to count the combinatorial 
factor explicitely, we need to prepare some notions concerning trees. 

Take any T G T(N„ + i) and for any q G N n +i let d q (g N) denote the incidence 
number of the vertex q (the number of lines connected to the vertex q) . From now 
let us always think that any tree in T(N n +i) starts from the vertex 0. For any 
q G N„+i let L q {T) (c T) be the set of lines from the vertex q to the vertices of the 
later generation. We see that 

tfL (T) = d , $L q (T) = d q - 1, Vq G N„ +1 \{0}. 

We define the combinatorial factor N(T) we want to calculate as follows. 

Definition 4.6. For any T G T(N n+ i) the combinatorial factor N(T)(e N) is 
defined as the total number of monomials appearing in the expansion of 

n (A g ,g/ + A q ,, q ) n ^,* Ta .j^ M « lx «^T*!^»ii*j- ( 4 - 13 ) 

Note that N(T) is independent of how to choose z q . e T, G [0, /?)/, (Vj G 
{1,2,3,4},V<ZGN„ +1 ). 

The combinatorial factor N(T) is counted as follows. 

Lemma 4.7. For T G T(N n +i) let d q {q G N n +i) denote the incidence number of 
the vertex q in T . If there is q G N„+i such that d q > 4, N(T) = 0. Otherwise, 

N(T) = i n ( d 3 

Proo/. Set := n ge N„ +1 V , z ?T^?^J|a|V , z |Tx^<;xl- If there is p G N n+i such that 
dp > 4, every term in the expansion of the product H{ q , q '}£T(A q , q ' + A q i jq ) contains 
more than 4 derivatives with respect to the Grassmann algebras indexed by p. Since 
the number of the Grassmann algebras with index p in W is 4, (|4. 13|) must vanish. 

Let us consider the case that d q G {1,2,3,4} for all q G N„+i. The operator 
A q . q i can be decomposed as A q ^ = X)cre{t n ^q,q'' where 

A ^ : =-E E C h ( X *x,y*y) 9 q * (4.14) 
x,yer x ,ye[O,0) h °^ ax Wyvy 

for a G {T, 4}- Note that for any a G {f, 1} and p,p',p" G N n+ i 

A^A^W = A;, jP A;„ iP ^ = 0. (4.15) 

By changing the numbering of vertices if necessary we may assume the following 
condition on T without losing generality. 
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(Jit) The distance between the vertex p and the initial vertex is less than equal to 
that between the vertex q and the vertex for all p,q € N„+i\{0} with p < q. 

Note that 

n (A q , q ,+A q ,, q )w= n n (a^+a^+a^+a^jw. (4.16) 

Let us count N(T) recursively with respect to q E N n +i as follows. The expansion of 
the product U{o, p }&l (t) ) is a sum of 4 ttL °' T ) terms, each of 

which is a product of §Lq(T) Laplacians. By the property (|4.15p any term containing 
the products A^A^, or AJ i0 AJ j0 for some a G {0, g}, {0, g'} e L (T) 

does not contribute to the number of remaining monomials in (|4.16[) . thus, can be 
eliminated. Therefore, we only need to count 



4 



ttio(T)! 



terms in the expansion of n { o,p}eL (T) ( A l, p + &o, P + A jLo + A p,o)- 

Take any q S N„+i\{0} with L q (T) ^ 0. By the condition (Jfc) there uniquely 
exists q' € N„+i with q' < q such that {(/',<?} € L q i(T). Thus, every term in the 
expansion of the product 

n n (A T + a 1 +A T - + A J --) 

3&f n+1 ,L j (T)^» {i, P }eZj(T) 3,v 3 ' p p ' 3 p ' 3 

contains one of the Laplacians of the form A q q , , A q , (cr € {f, J.}). Therefore, by 
the property (|4.15p only 

3 

$L q (T) 

terms in the expansion of the product H{q, P }eL q (T) (Aj )P + A^ p + A^ p q + A p ? ) need 
to be counted. 

By repeating this argument recursively for all q <G N n+ i\{0} we can calculate 
N(T) as follows. 

N{T) = ( ttXoV) ) » i0(T)! ( WT) ) «^ (T)! 



tti,(T)! 



4 



do! n , " (d,-l)! = 4 n " (d,-l)! 

<jeN„,,\(ov \ a n — l / ofM.x, \ a n — 1 / 



<^0 / geN„ +1 \{o} y ««j — 1 J ' ' geN„ + i \ u q 

d q >2 

where we used the fact that 



0! = 1. 



□ 
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Let Dh denote the L^norm of the covariance matrix Ch, i-e, 



Dh: =jY. E |Ch(xTa: s 0T0)|. 



x£riE[-W) h 

Now we can find an upper bound on |oh in |. 

Lemma 4.8. For any h S 2N//3 and for all n e N U {0} i/ie following equality 
holds. 



\ah,n\ < 



128 

3n + 4 



3n + 4 
n 



(41> h ) n . 



(4.17) 



Proof. By using Lemma [2T5l and (|3.2[> - (|3.3p we can directly evaluate \a,h,o\ to obtain 
\ah,o\ < 32, which is (fCT7|) for n = 0. Let us show (|4T7|) for n > 1. By flU]) and 
(|4.8|) . we need to evaluate 



1 

/to! 



V n (A q , q ,+A q ,, q )[ ds V 0(T,7t,b) 
T(N„ +1 )^' 9 '> eT -Wl" .es„ +l( T) 



,A(Af(T,7r,s)) 



z„er 



(4.18) 



By using the property l|4.9p we have 



E 



j I / sup 

pn! f-^ se[o,i]" 



3 A(Af(T,7T,s)) 



r n (A, ja /+A,/, 9 ) 

{?,9'}eT 



•*Wi&,i&) n E*W^,^) 



z„er 



(4.19) 



Take any T 6 T(N„ + i). If T contains a vertex whose incidence number is larger 
than 4, by Lemma T4. 71 the tree T does not contribute to the sum in (I4.19|) . Thus, 
we may assume that the incidence numbers do, d\, ■ ■ ■ ,d n of T are less than equal 
to 4. Moreover, as in the proof of Lemma 14.71 without losing generality we may 
assume the condition (Jfc) on T. 

Let qi, <72, ■ ■ ■ ,9; G N n +i\{0} be such that qi < q2 < ■ ■ ■ < qi and 
{<?i> 92, • ' ' , <li} = {<? S N n+ i\{0} | L q (T) ^ 0}. Every term of the expansion of the 
product ri{ g , g '} G T(^9,9' + has the form 



n a? -, n n a? ,, 



(4.20) 



where by using the notation (HTH . A^ p} € {A T qp , A T pq \ r € {T, I}} for all e 
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L q {T) with q £ N n+ i satisfying L q (T) ^ 0. For any {0,p} £ L {T) we set 

C{o,p}( x O' z p x p) '■= 

|Cfc(z p Ta;p,xita;o)| if A {o, P } = A P ,0' \Ch(z P lx p , x 2 |x )| if A{ 0p} = A p0 , 
|C^(yitxo,ZpTa:p)l if A {o :P } = A o,p' \Ch(y2ix ,z p ix p )\ if A^ p} = A^ p 

for all xq, x p G [0, 0)h and z p 6 T. For any j G {1, 2, • • • , 1} and any {qj 7 p} G i 9j (T) 
we define 

C {qj,p}(?<i3 x <iii*p x p) : =\ C h(z qj TX qj ,z v TX p )\ if A{^, p} = A^ iP for some r G {T, I}, 
|C h (z p rXp,z gj Ta; g3 )| if A^. p} = A£ g . for some t g {f, j} 

for all x p ,x qj £ [0,/3)/j and z p ,z 9j G I\ By noting that (14.20|) is the product of n 
Laplacians and using Lemma [4751 and the condition (Jfr), we observe that 

e A ( M(T,., s) ) n A ? 0po} A n A? , 



v h ,zM°x,i?x) n E v^Mxrfx) 

9=1 — L 



< 



z„er 



E n (~ E e) n 



C{o, Po }( x °> Z Po x Po) 



l> ^ {o P oKMT) U S E Gfo, P „}(^o,z Po .T Po ) 



A71 + 2 



{9i,pi}ei 



{?2,P2}6i 



{91 ,pi}ei 



n I T E E C {quPl}( Z 1l X q^ Z Pl X Pl) J 

5 r ^(7; E E C 'te.P2}( Z 92 :Z; «2 J Z P2 :Z; P2) ] 

ei ^ (T) V x P2 e[o,p) h * P2 er J 

n (i E E 



( z 9i 2-9! 1 Z P( x Pi ) 
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n Sup - ^2 C {0,p }( X 0> Z P« X Po)\ 

'- z « er Ipi e[o,/3) tZpi er / 

T J]] C{g 2 ,p 2 }( Z g2 2; ?2J Z P2 :Z; P2) J 

x P2 e[o,/3) h z P2 er / 

SUP ^ X! X! C {qi,pi}( Z H X 1l' Z Pl X Pl) ) 

« l e[o ) « h ,«„erA !i!pi6[0i/S j hBpj6r y 



< 4 n+2 /3 

{0,po 



n sup 

{«i,pi}ei 91 (T) U w e[o,/3)h, 



n sup 

{92,P2}eL, 2 (T) \ £c„ 2 e[o,/3) h ,z 92 er 



n 

{?i,Pi}ei 9! (T) 
< 4 n+2 /?£>£. 

(4.21) 

To obtain the last inequality in (|4.2ip we have used the fact that for all {0,po} € 
L (T) and all {q 3 ,p 3 } G L qj (T) (Vj e {1, 2, • • • , I}) 



SUP \ Z Z C {0,p }( x O,Zpo x Po) < D h, 



By combining (|4. 2 1 [1 with (|4. 19|) we have 



i « 

1S3<4W E N(T), (4.22) 

T£T(N„ + i) 



where N(T) is defined in Definition 14.61 

By repeating the same argument as above we can show the inequality (|4.22p for 
the case that Zq — X%, Thus, by recalling (14. ip we arrive at 

\a h , n \<-M^hT V N(T). (4.23) 
n! ' 

TeT(N„+i) 

To complete the proof we calculate the sum Stst(n„ + i) N(T) in (|4.23[> . As 
characterized in Lemma |4.7[ the number N(T) only depends on the incidence num- 
bers of T. By using Lemma 14.71 and Cayley's theorem on the number of trees with 
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fixed incidence numbers (see, e.g, (22 Corollary 2.2.4]) we have 



E N(T) 

TeT(N„ + i) 



E 1 E, ! 6K Il+1 d q =2n 



(n-1)! 



d,6{l,2,3,4} 
Vg£N„ +1 



- • 4 n 

U q en n+1 (d q - 1)! <?eN„ + i 



K - 1)! 



4(n-l)! V l v d 



d,e{l,2,3,4} 
V 9 EN„ + 1 



.an- n 



Moreover, by using Cauchy's integral formula and the residue theorem we see that 
for a positive r > 



E N(T) = 

TGT(N„+i) 

_ 4(n- 1)! 



4(n- 1)! / d 



(2n)l 



2)1 



E 

v.d=l 



3 

d- 1 



n+l 



2?ri 
4(n- 1)! 



|js|=r 

3ra + 3 
n — 1 



z n+l/1 i ,\3(n+l) 

dz = 4 (« - 1)! Res 

z 2n+l V ' Z=Q 



z=0 

(1 + z) 3 (" +1 ) 



3n + 4 



3n + 4 
n 



Combining (|4~24|) with (|4~23]) yields the result. 



(4.24) 
□ 



The inequality (14. 17)) motivates us to know the properties of the power series 
f(x) defined by 



/(*) ■■= E 



3n + 4 



3n + 4 
n 



(4.25) 



As the last lemma before our main theorem, the properties of f(x) are summarized. 

Lemma 4.9. The radius of convergence of the power series f{x) defined in (|4.25[) 
is 4/27 and /(4/27) = 81/16. Moreover, for any x € (0,4/27] the following equality 
holds. 

(* 



16 

9x2 



tan 



-± 1 

27x 



(4.26) 



V 



where the function tan 1 (-) is defined as a bijective map from R to (— 7r/2,7r/2) 
satisfying tan _1 (tan(6>)) = 8 for all 9 E (— w/2, 7r/2). 

Proof. As a topic in generating functions the power series (|4.25p is commonly stud- 
ied (see, e.g, pp. 200-201]). However, we give a proof for the statements for 
completeness. Let us analyze the cubic equation 



X = 1 + xX 6 



(4.27) 
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forxg (0,4/27). We see that for any x € (0, 4/27) and z £ C satisfying | z- 1| = 1/2, 
the inequality \xz 3 \ < \z—l\ holds. Thus, the Lagrange inversion theorem (see, e.g, 
[TBI Theorem 2.3.1]) implies that in the domain {z £ C | \z — 1| < 1/2} there is 
exactly one root X = v(x) of (|4.27[) and 



v(xy = i + j2 



El (A. 

n\ \ dz 



(4zV*) 



fix). 



(4.28) 



On the other hand, by algebraically solving (|4.27|) and specifying a root con- 
tained in the domain {z S C | \z — 1| < 1/2} we can determine the explicite form 
of v(x) as follows. 



v(x) 



: COS 



tan 



- 1 (v / S~ T ) + ^ 



(4.29) 



V 



where tan 1 (-) is defined as stated in Lemma T4. 91 

The equalities P~28]l - (|3~29"ll give ([426]) for x G (0,4/27). The ratio test shows 
that the radius of convergence of / is 4/27. Moreover, by continuity we have 



lim 



xS4/27 



f(x) = u(4/27) 4 = 81/16, which completes the proof. 



Define a constant D > by 
D := lim 



D h = I da ! y;|C , (xTa! > OT0)|. 
-P xer 



□ 



(4.30) 



Our main result is stated as follows. 

Theorem 4.10. For any xi,X2,yi,y2 € T and mSNU {0} and any U £ R with 
\U\ < 1/(2713), the following equality and inequalities hold. 



KdrCa^i^yiT + V' yi tV'y 2 iV'x 2 iV'x 1 T)l < ^(It 7 !), 

m 
n=0 

128 / 3n + 4 



< 



R(\u\)-J2 



3n + 4 



(4.31) 
(4.32) 

(4.33) 



(AD\U\) n , 



where {a n }^L Q is given in (|2.11[) - (|2.12[) and 

32 

tan - 



R(\U\) := 



32 



27D\U\ 



- 1 



9L> 2 |[/| 2 



ifU = Q, 
if0<\U\<^, 



(4.34) 
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with the function tan : M — > (— ir/2, tt/2) satisfying tan x (tan0) = 9 for all 

e e (-7r/2,7r/2)- 

Proof. Since by Lemma 13.11 (fm|) and Lemma 14.81 



K| = lim |o M | < ( 3n + 4 ) (47?)", (4.35) 

Lemma SU implies that for all [7 G [-1/(271?), 1/(271?)] 

OO 

J2K\\U\ n < R(\U\), (4.36) 

n=0 

where i?(|C/|) is defined in lpL"5g]) . The inequalities l[02]) - l[g^5)l follow from (|OT|) 
and (|L35 l) - (jt36l) . 

We show that the equality P~3"l"]) holds for U G R with |C/| < 1/(27|7?|). Let us 
fix any e G (0,1/(271?)). Since P\x x=0 ,vxer* = Tr e"^/ Tr e"^ > for all [7 G 
K, Proposition 13.21 implies that there exists Nq G N such that \Ph\x x =o.\/xer i \ > 
for all h G 2N//3 with /i > N //3 and all [7 G K with \U\ < 1/(271?) - £. Moreover, 
since Ph\\ x =oyx£T 4 is a polynomial of U we can take a simply connected domain 
Oh (c C) containing the interval [—1/(271?) + e, 1/(271?) — e] inside such that 
\Ph\\ x =oyxerA > for all U G O h . Thus, we see that dPn/dX^ / Ph\\ x =oype&T l 
defines an analytic function of U in the domain Oh- By Lemma 14.81 and Lemma 14.91 
the series Y^=o a h,nU n converges for all U G C with \U\ < 1/(277?^). By choosing 
7V sufficiently large we may assume that 1/(271?) - e < 1/(277?,,) for all h G 2N//3 
with h > No/ (3. Therefore, Lemma [3.11 JO]) and the identity theorem for analytic 
functions ensure that 



1 dT~ P h({ u x}xer^) 



\ x =o 



= J2a h , n U n (4.37) 



n=0 



Ph{{U x }x^) 

for all [/ G [-1/(271?) + e, 1/(277?) - e]. 
Note that Lemma 14.81 implies 

rrn . 128 /3n + 4\/4\" 
I^ ^^Tii n J (,27 J (4 ' 38) 

for all U G [-l/(277?)+£, 1/(271?) -e] and the right hand side of is summable 

over N U {0}. Thus, by Lemma f3.1l IpS]). Corollary 13.41 and Lebesgue's dominated 
convergence theorem for / 1 (N U {0}) we can pass h — > +oo in (|4.37p to deduce the 
equality IpTBT) for all [7 G [-1/(277?) + £,1/(277?) - e]. Then, by sending e \ 
and continuity we obtain (|QT|) for all [7 G [-1/(277?), 1/(277?)]. □ 



Remark 4.11. In Proposition 1 5 . 1 1 we give a volume- independent upper bound on 
the decay constant 7? in 2 dimensional case. One can straightforwardly extend the 
calculation of Proposition 15.11 to derive volume-independent upper bounds on 7? 
in any dimension. By replacing 7? by these upper bounds, Theorem 14.101 provides 
volume-independent upper bounds on the perturbation series X^i^Lo a nU n ■ 
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5 Numerical results in 2D 



In this section we compute the perturbation series of the correlation function 
(V^t^i^WVViT + V'y lT '0y 2 |'0x 2 i'0x 1 T> UP to 2nd order ter m in 2 dimensional 
case. We also implement the upper bound obtained in Theorem 14. 101 and report the 
error between the correlation function and the 2nd order perturbation. Throughout 
this section it is assumed that d = 2. 



5.1 The decay constant for d = 2 

In order to estimate the radius of convergence of the perturbation series E^Lo a nU n 
and compute the upper bound on the sum of the higher order terms numerically, 
first we need to evaluate the decay constant D defined in (|4.30|) . The result is 
presented in the following proposition. 



Proposition 5.1. The following inequality holds 
D < 



(5.1) 



16 32tt 2 16tt 3 \ lfP_ (2jt|+4jt])e^g 
W + 3Vfp + ) \ 2 + i + rfi 

o ( (2\t\+mv^ \ 2 , f m+m^i x 

+ 6 \ 1 + eK J + { l + e K ) 

where £ := 4|t| + 4|i'| + 

The derivation of the inequality (|5.ip needs the following estimate. 

Lemma 5.2. The following inequality holds. 

1 , 8tt 2 P 4tt 3 /3 2 

> < 4 H — H — 

Ztr 1 + E?=i \<? 2 ™ l/L - 1| 3 L 3 /(8tt 3 /? 3 ) " 3^3 3^3 

where x = (2:1,2:2) G T. 

Proof. For any y £ R let |_2/J denote the largest integer which does not exceed y. 
By using the inequality that \e l6 — l| > 2|0|/7r for any 6* 6 [— 7r, 7t], we see that 

y I 

Z^r 1 + Ef=i le^'A - l| 3 L 3 /(8^ 3 /? 3 ) 

LL/2J LL/2J 



4 E E T , ^ 



Xi= 

OO OO 

< 



4 E E i 



o 1 + Ej=i |e i27ra: '/ L - 1| 3 L 3 /(8tt 3 /3 3 ) 
1 



„ „ - + 8x 3 /(7r 3 /3 3 ) + 8a; 3 /(7r 3 /3 3 ) 



i ^ oo oo 

E l + 8x 3 /(7r 3 /3 3 ) + 4 E E i 



+ 8x 3 /(7r 3 /3 3 ) + 8x 3 /(tt 3 /3 3 ) 
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4 " 8 /o ^l + S,?/^) 



poo poo 

4 / dxi / dx2 



1 + 8x 3 /(tt 3 /3 3 ) + 8x 3 /(tt 3 /3 3 ) 



= 4 + 4-7T- > / das - - ^ + tt 2 /3 2 / <l.r 



(l + X 3 )2/3 7 l + a; 3 



oo j 

dx- 



< 4 + (4ir/3 + 2tt 2 /? 2 ) / da:- ?> 

Jo 1 + ^ 



1 



where we used the inequality 

/■OO 

J ^'(l+X^ 

Then, by using the equality 



< 1 



= 2. 



/ dx- 

Jo 1 + - 



2tt 
3^3 



(see [TU]), we obtain the desired inequality. 



□ 



Proof of Proposition I5.il Let us define a linear operator di l ■ C°° (R 2 ) — * C°° (R 2 ) 
(1 = 1,2) by 

/(k+27rei/L)-/(k) 



(<W)(k) := 



2tt/L 



for any / 6 C°°(R 2 ). Then, the mean value theorem shows that for all k = (fci, k 2 ) £ 
ra>2 



|(di,L) 3 /(fci,fe)| < sup 

fciS[fci,fci+67r/L] 

|(d 2 ,L) 3 /(fci,fe)| < sup 

*:2e[*;2,fc2+67r/L] 



a 3 - 

<9 3 



(5.2) 



Define a function F(k, x) : R 2 x R -> R by F(k, a) := e^/C 1 + e /3£ k ). We see 
that for any x = (x\, x%) £ T, x £ [—0, 0] and I £ {1,2} 



2-k/L 



C(xtx,OT0) 



l_ ^2 e-^((d liL ) 3 F(k,x)l x > - (d liL ) 3 F(k,x + I3)l x<0 ). 



(5.3) 



I'' 



ker* 



By (fS~2"l) and (|5~5|) we have that for any x = {xi,x 2 ) £ T, x £ [-(3, (3] and I £ {1, 2} 



,i2irxi/L _ 

2n/L 



C(xtx,OT0) 



< sup sup 

se[0,(9] ke[0,2ir+67r/L]x [0,27r+67r/L] 



a 3 
9fcf 



F(k,a;) 



(5.4) 
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Note that for I = 1,2 
d 3 , d 3 E k d 



„xE k 



dE k d 2 E k d 2 



(dE k \ 


3 Q3 


( e xEk \ 


\ dh ) 




V 1 + eP E * ) 



(5.5) 



and 



|£ k |<4|t|+4|t'| + | M |, 



dE k 



dk. 



d 2 E k 



dk 2 



d 3 E k 



dk 3 



<2|i|+4|i'|. (5.6) 



Moreover, we can prove that for m S {1, 2, 3} and any E G 



sup 

z€[0,/3] 



d 

dE 



xE 



< 



l + e^l J 



(5.7) 



By the equalities that e xE /{l + e pE ) = e^- x )(- E )/(l + e^- E )) and d/dE = 
—d/d(—E) we may assume E > without losing generality in the following ar- 



ument. First 


let us 


study 


the 


case for m 




d 






e xE / 




dE 1 


+ e 0E 




1 + 



0e 



.dE 



l + eP E 



< 



0e 



0E 



l + e^' 



which is (|5.7|) for to = 1. 

Next consider the case that m — 2. We see that 



„xE 



^xE 



-9i{x): 



with 



dE J l + eP E l + e^' 



gi(x) := a; 



l + e^ 



(l + e^) 2 



Note that 101(0)1, \g\(peP E /{\ + e^ E ))\, \gi(ff)\ < 2 e^ E /(l + e' iE ). Thus, we have 
that 



d 2 e XoE 



dE 2 1 + eP E 



< 



< 



l + e? E 



max{M0)|, \ 9l (^ E /(l + e^))|, | 5l (/3)|} 



I3E-, 



Be 



l + e0 E , 

Finally, analyze the case for m — 3. A calculation shows that 



xE 



a xE 



dE l + eP E 1 + e? E 



92(x), 



(5.8) 



where 

g 2 (x) := x 3 



0E 



3/3e 



3/? 2 (e^-l), 

(l + e /3B)2 



(l + e^) 3 



37 



The roots x\, x 2 of the equation g' 2 {x) — are given by 

Since < x\ < (3 < x 2 , max ie[0i/ j] \g%{x)\ = max{|g 2 (0)|, \g2{%i)\, \92(P)\}- More- 
over, we see that 

Q3 p /3E ftZpPE fl3 p 2/3E 



( 1 + e /3£)3' ' - (1 + e l3E)3 V / ( 1 + e /3B)2' 

l32(/?)l = (1 / e 3 ^ )3 l e2/9£ " 4e/5B + !| ^ M°)l> 

fl3„/3B f)3 p /3E ft3 p 2f3E 



pi2-KXi/L _ ^ x " 

2tt/L 



C(xtx,OT0) 



< (2|t|+4|t'|)/?e^ | 3 /(2ltl+4lf'l)/?e^V , ( (2|t| + 4|f|) J 9e« 



(l + e^) 3 ' 1 " (1 + e^) 3 v ' ' (l + e^) 2 ' 
which imply that 

p3 e 20E 

^ Ma0)| - (1 + ^ - (5J) 

Combining (JSTHJ) with |Q| deduces (|57]) for m = 3. 
Then by ([Q] ) -([5?f |l we have for Z = 1, 2 



(5.10) 

1 + eK 1 " ^ 1 + e /3« ; ' V 1 + eK ) ' 

with £ := 4|t| + 4|i'| + By using the inequality that |C(xta;, 0|0)| < 1, IpTTOjl 
and Lemma 15.21 we can derive the inequality (|5.1[) . □ 

5.2 The 2nd order perturbation 

As a preparation for the numerical implementation of the 2nd order perturbation 
ao + a\U + a 2 U 2 , we rewrite the formula (|2.11[) for a n in a more suitable form for 
practical computation. 

By decomposing the determinant of the covariance matrix into the sums over 
permutations, G n defined in (|2.12p can be written as follows. 

G n = — V sgn(7r)sgn(T) 5 „(7r,T), (5.11) 
n! 

W,T£Sn 
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with 



n 



xi,x 3 2 ,y 3 1)y ^er 







n"xi,x 3 2 V 1 ,yrxi,y 3 I ^ 'M.x'.yi.y* T" - \v .y ,.x .x, 



C(xj Txj , y[ w Tavy) )C(x 3 2 J.xj , yJ w K(j) ) 



The formula (|2.11[) for ao,a\,a2 is 
1 d 



a = 



«2 



i a 



Gi 



i d ( 1 2 



G3 — G1G2 + — 



a a .=o ' 



(5.12) 



By substituting (|5.11[) into (|5.12p and canceling we can simplify the expressions for 
a\, a,2 as follows. 



«i 



19 1/ 



1 2 

2 1 



1 2 

2 1 



ft Id 



« 2 



+ .93 

- 353 

+ .93 



1 2 3 

2 3 1 

12 3 
1 3 2 

1 2 3 

2 3 1 



where 



Id 2 
M3 



I 




i 






n 




n 



1 2 3 
1 2 3 



1 2 

2 1 



1 2 3 

2 3 1 

1 2 3 

2 3 1 

12 3 
2 3 1 

1 2 3 

3 1 2 



G S 2 

1 2 3 
1 3 2 



.92 

"53 

h3 33 
- 3^3 



A^=0 

v^er 4 



1 2 

2 1 



,M ; 



x x =o 



1 2 3 

2 3 1 



,Id ; 



: 


2 




i 


2 






3 






1 





1 2 3 

2 3 1 



1 2 3 
1 3 2 



1 2 3 

2 3 1 



1 2 3 
3 1 2 



G S 3 



In practice we need to implement the permutation-dependent terms 
d/dX^ g n (ir, t)| A^=o.VA'er 4 - Let us describe its implementation below. 
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Note that the multiple integral n" =1 Jjf dxjf{x\, • • • , x n ) of any integrable func- 
tion / can be decomposed into a sum of nl integrals as follows. 

n [P n [P 

II / dXjf(xi,---,X n )= II / d,Xj 2, li,(i)>x,(2)>'->^(„) /( s l i ' • ' i ^n) 

3=1 Jo i=Uo neSn 



En f x V<J-V 
II / dx v <j)f(xi, ■■■ ,x n ), 

a 3=1 Jo 

(5.13) 



(5.14) 



where £E„(o) := /3. 

By substituting the expression (|2.7[) of the covariance matrix and decomposing 

the integral fl™ =1 J dxj as in (|5.13[) , we can expand d/dX^ i g n (n, T)| AA , =0 ,v;ter 4 as 
sums over the momentum space T* as follows. For neN 

g , , _ (-1)" \- 

f)\. 9n\T^,T) Xx=Q - {n+1)d 

* v * er4 k 1 ,...,k n , P i,...,p B er' 

n 

■ 2j(cos((k,-,xi) + (pj,x 2 ) - (k^-i^yi) - (p T -i (j - ) ,y 2 )) 

3=1 

+ cos((kj,yi) + (Pj,y 2 ) - (k^-iy), x x ) - (p T -i (j) ,x 2 ))) 

n 

• n <5(k z + p; - k^-i^) - p T -i(j)) 

i?y 

r,esJ= lJo 

n f lx 7 , u) -x J <a _ la: j(j) -Xj>o \ / lay^-a^O _ ^x TU) -Xj>0 

' /il ^ i + e ^ kj i + e ~0E kj + e 0E Pj 1 + e -0E Vj 

where x„(o) — 

Let us sketch how to implement (|5.14p . We prepare a function of real vari- 
ables Ei, ■ ■ ■ ,E n returning the exact value of II™ =1 Jq 31 dxje XjEj with xq = 
beforehand. Then we iterate the system with respect to the variables ki, • • • , k„, 
Pi) • • " > Pn € T*. For fixed ki, • • • , k„, Pi, • • • , p n <E T* we iterate with respect to 
the permutation r/ G S n . For each r\ € S n we substitute the variables Ej = E^ ^. + 

E„ ... — Eh , — E , (j = 1, n) into the function IT?-i f„ Jl dx 
and its returning value is then multiplied by the constant 



p X 3 

3 C 



^-x^t^ —Xj<0 ^-x 7T (j)—Xj>o\ ( ^-x T tj) — Xj <0 lavm — xj >0 



3=1 \ 1 + e ■> l + e j/\l + e^ p J 1 + e p 
5.3 Numerical values 

Here we display our numerical results. In our computation we fix the physical 
parameters t,t',fi,/3 to satisfy t = t' = ji = 0.01, j3 = 1. In this configuration the 
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\u\ 


1.0 x 10~ y 


5.0 x 10"° 


1.0 x 10~ 5 


5.0 x 10~ 5 


Error 


1.408 x 10~ 7 


1.773 x 10~ b 


1.433 x 10~ 4 


1.942 x 10~ 2 


\u\ 


1.0 x 10~ 4 


2.0 x 10~ 4 


3.0 x 10~ 4 


4.0 x 10~ 4 


Error 


1.739 x 10- 1 


1.842 


9.454 


7.307 x 10 



Table 1: Errors between the correlation function and the 2nd order perturbation 



L 


10, ll,-- - ,18 


a 


5.050 x lO" 1 


Ol 


-3.774 x lO" 1 


a 2 


9.339 x 10- A 


a + aill + a 2 U 2 


5.050 x 10" 1 



Table 2: 2nd order perturbation in the case that Xi = x 2 = yi = Y2 = (M) 

(le{o,i,---,5}) 



upper bound on D obtained in Proposition 15. II is 92.04. Thus, by Theorem 14. 1 01 the 
radius of convergence 1/(27|D|) of our perturbation series X^o a nU n is estimated 
to be larger than equal to 4.024 x 10~ 4 . 

The errors between the correlation function and the 2nd order perturbation for 
various \U\ less than 4.024 x 10~ 4 are exhibited in Table [TJ where Error is defined 
by the right hand side of the inequality (14. 33[) for m — 2 and D = 92.04 and satisfies 
that 

KdrCa^-lAiT + V'y lT V'y 2 iV'x 2 iV'x 1 T) - ao - aiU - a 2 U 2 \ < Error. 

Let us fix U = l.Ox 10" 5 . According to Tabled] the error between the correlation 
function and the 2nd order perturbation is estimated as 1.433 x 10 -4 . Table[5]shows 
values of a , a\, a% and a + a\U + a 2 U 2 in the case that Xi = x 2 = yi = y 2 — {I, I) 
(I G {0, 1, • • • ,5}) for various lattice size L from 10 up to 18. We observe that 
each of ao, ai, a 2 respectively takes the same value for any L € {10, • • • , 18} and 

le{0,--- ,5}. 

Table [3] shows the values of ao, Oi, a 2 and ao + a%U + a 2 U 2 in the case that 
xi = yi = (0, 0), x 2 = y 2 = (1,1) (I E {I, ■ ■ ■ , 5}) for various lattice size L from 10 
to 18. Again we see that each of ao, a\, a 2 respectively takes the same value for any 
L e {10, • • • , 18} and I € {1, • • • , 5}. Since we have fixed a small U so that Error 
becomes sufficiently small, the 1st and 2nd order terms do not contribute to the 
sum ao + aiiJ + a 2 U 2 much in these numerical simulations. 

We also computed 00,01,02 in the case that xi = x 2 — (0,0), yi — y 2 = (1,1) 
for I e {1, • • • , 5} for L = 10, 11, • • • , 18. The result shows that |oo|, |oi|, 1 0,2 1 < 
1.5 x 10 -5 for any i e {1, • • • ,5} and L G {10, 11, • • • , 18}. In this case the values 
of |oo|, |oi|, I02I are much smaller than those presented in Tables [2][3] This result 
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L 


10,11,- ■■ ,18 


a 


5.050 x 10" 1 


Ol 


-2.524 x lO" 1 


a-2 


9.402 x 10~ 2 


a + aiU + a 2 U 2 


5.050 x 10- 1 



Table 3: 2nd order perturbation in the case that Xi = yi = (0,0), = y? = (1,1) 
{I €{!,•■■ ,5}) 



indicates that the 4 point correlation function takes small values if xi = X2 , yi = y2 
and |xi — yi | is large and agrees with the decaying property of the 4 point correlation 
function for the 2 dimensional Hubbard model proved in 12J. 



Appendices 

In this section we review the definitions of the Fermionic Fock space and the an- 
nihilation, creation operators, prove Proposition 12.41 and show that the covariancc 
matrix Ch has a non-zero determinant independent of the parameter h. 

We write a matrix M indexed by finite sets S, S' with §S — jjS" = n as M = 
(M(s, s')) S £S.s'eS' ■ I n this notation let us think that each element of S and S' has al- 
ready been given a number from 1 to n and M is defined by M — (M(sj, s' k ))i<j t k< n 
even if the numbering of S and S' is not specified in the context. The main re- 
sults Proposition 12 .41 and Proposition IC . 71 concluded after some argument involving 
such matrices in this Appendices are independent of how to number the index sets. 
For any finite set B let L 2 (B\C) denote the complex linear space consisting of 
complex- valued functions on B, even when we do not introduce an inner product 
in L 2 (B; C). 

A The Fermionic Fock space 

The first part of Appendices reviews the definitions of the Fermionic Fock space on 
the lattice T x {t, 1} and the annihilation, creation operators ^xj,^. 

For any n £ N we consider the linear space L 2 ((T x {T,|}) n ;C) as a Hilbert 
space equipped with the inner product (•, ')L 2 ((rx{T 1})™ C) defined by 

(01) <h) £2((rx{ta}) n ;C) 

:= <Mxi<7l, • • • ,X„CT„)0 2 (xiCTl, • • • ,X„<7„). 

xi,--- ,x„er CTl ,--. ,<r ra e{T,|} 
By convention we set L 2 ((r x {|, |})°; C) := C. 
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For n E N the anti-symmetrization operator A n : L 2 ((T x {|, |}) n ; C) — » 
i 2 ((r x {t, J.})"; C) is defined by 

(A„0)(xicti, • • • ,x„ct„) ^ sgn(7r)0(x 7r(1)< T 7r ( 1 ),- • • .x.^^^). 

The operator is defined as the identity map on C, i.e, Aqz := z for all z € C. 

The subspace A„(£ 2 ((T x {!, |})";C)) of L 2 ((T x {T,|})";C) is called as the 
Fermionic n~ particle space and is a Hilbert space equipped with the inner product 
of L 2 ((T x {T, |})™; C). Note that by anti-symmetry for any n > 2L d , 

A„(L 2 ((rx{t,|})";C))-{0}. 

The Fermionic Fock space Ff(L 2 (T x {|, J,}; C)) is defined as the direct sum of 
4(i 2 ((rx{t,I})";C)) (n = (),••• ,2L d ) as follows. 

2L d 

F f (L 2 (T x {T, |}; C)) := A n (L 2 ((T x {|, {})"; C)). 
The space Ff(L 2 (T x{|,J.};C))isa Hilbert space with inner product (•, -) Ff defined 

by 

2L d 

(01, fa) Ff ■= ^ ^ 1 >™' ^2.™)L2((rx{T, !})";€)' 

Tl = 

for any vectors fa = {fa,o,fa,i, ■ ■ ■ ,fa,2L<t), fa = (fa.o, fa,i,- ■ ■ ,fa,2L d ) 

e^(L 2 (rx{T,|};C)). 

Define a set of functions {faa}(k,a)er*x{t,i} C i 2 (r x {T,|};C) by 
<^k<r(xr) := S a T e~ %1 ^-^ jL d l 2 . We then define a function 0ki<n ■ ■ ■ fa n a n 
GL 2 ((rx{t,|})";C)by 

0ki CTl ' ' ' 0k„<r„(xiTi, • • • ,X„T„) := 0k 1<Tl (xiTl) • • ik„a„(x n T„). 

An orthonormal basis of Ff(L 2 (T x {|, J.}; C)) is given by Un=o wnere 
So :={1}(C C), 

B n := < V^!A„ [ n fa^J j n kCT G {0, 1}, V V n k(J = n I 
I Wr*x {U} y ^,em } J 

(A.l) 

for n e {1,2, ■ ■ ■ , 2i d }. Thus, we see that dimf>(£ 2 (r x {|, |}; C)) = £^1 j}5„ = 

The annihilation operator Vx<x : F f {L 2 (T x {t,j};C)) -» F f (L 2 (T x {t,j};C)) 
(x e r, a € {t, I}) is defined in the following steps. For any n G N U {0} and any 
4>eA n+1 (L 2 ((rx {U})" +1 ;C)), ^<f> e A n (L 2 ((r x {U})";C)) is defined by 

(V'xo-^CxiCTi,- • • ,x„cr„) := 10(xcr,xicri,- • • ,x„cr„). 
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For any z £ Aq(L 2 ((T x {f, J.}) ; C)), "0xo-2 := 0. The domain of the operator i/) xa 
is then extended to the whole space Ff(L 2 (T x {"f , J,}; C)) by linearity. 

The creation operator : Ff(L 2 (T x {t, |}; C)) -> F f (L 2 (T x {|, J.}; C)) is the 
adjoint operator of -0xcr and characterized as follows. For any ngNU {0} and any 

0e A„(L 2 ((rx {U})";C)), ^> e A„ +1 (i 2 ((r x {U})" +1 ;C)) and 

(V'x O -0)( x i cr ir • • ,x„ + icr„ + i) 
, n+l 

= , V(-l)' _1 '5x,x i '5 CT . (T! (/)(xicri, ■ • ■ ,5q07, • • • ,x n+ icr„ + i), 

where the notation 'xTct;' stands for the omission of the variable x;<7/. 

For any operators A, B let {A,B} denote AB + BA. The operators ipxa, Vxo- 
(x G r, a £ satisfy the following canonical anti-commutation relations. For 

all xjer, f7,re{T,l}, 

{V'xa, V'yr} = {V&rl lpy T } = 0, {V'xrr, V'yr} = <*x,y£<r,T- ( A ' 2 ) 

See, e.g, [3] for more detailed definitions of the Fermionic Fock space and the oper- 
ators on it. 



B The temperature-ordered perturbation series 

In this section we present the derivation of the perturbation series (|2.6p . Proposi- 
tions claimed here are standard tools in many-body theory (see, e.g, 
[TBI Chapter 2, Chapter 3]). This part of Appendices is devoted to show them in a 
mathematical context. 

Let us fix notations used in the analysis below. Let Hq, V\ be the operators de- 
fined in (|2.1|) - (|2.2p and (|2.3[) . respectively. In our argument in this section, however, 
we do not use the relation (|2.4p or the condition ()3.16p imposed on the parameter 
{Ux}xer 4 - One can consider more general V\ of the form (|2.3p parameterized by 
any complex multi-variable {Ux} a'gt 4 m this section. 

Define the operators V\(s), ipx^s), ^* CT (s) {s £ M, x £ T, a £ {T,4}) by 

V x (s) := e sH °V x e~ sH °, ^ a (s) := e sH ^ xa e- sH ° , i>* a (s) := e sfl »C/" Bl 
For a £ {0, 1} the operator ip xaa {s) is defined by 

if o = 0. 

Next we define the ordering operators T±, T%. 

Definition B.l. Consider linear operators M(si), • • • , M(s n ) 
: Ff {L 2 (r x {T,|};C)) -» F/(L 2 (r x {?,!};€)) parameterized by si,-- - ,s„ € 
K. Assume that Sj 7^ Sfc for any j, k £ {1, ••• ,n} with j 7^ fc. The operator 
Ti(M(si) • • • M(s n )) is defined by 

Ti(M(«i) • • • M («„)) := M(s ff(1) ) • • • M(« w(l0 ), 
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where ir € S n is uniquely determined by the condition that 

> s tt(2) > ' ' ' > S^( n ). 

Let us define a relation '~' in the set {V'x ( To(s)}(x, < 7,o, s )erx{T,i}x{o,i}xM as fol- 
lows. 

VWi(si) ~ ^yra 2 (S2) if Gtl = d 2 and Si = S 2 . 

We see that '~' is an equivalence relation in {Vw(s)}(x,o-,a,s)erx{T,.l}x{o,i}xiR- Let 
[V'x(to(s)] denote the equivalent class represented by an element Vxo-a(s)- We define 
relations V and 'h' in the quotient set {^o(s)}(x,<7,a, s )erx{T,I}x{o,i}xi/ ~ as 
follows. 

b/WiOi)] y bPyTa 2 (s2)] if si > s 2 , or Si = s 2 and oi > a 2 , 

[VWx^l)] h [V> yra 2 yra 2 (s 2 )] or h/> yra 2 (*2)]- 

The set {Vw(s)}(x, CT ,a, s )erx{T,.l}x{o,i}xR/ ~ is totally ordered under the relation 
'>r' and the relation V is a strict order in this quotient set. 

Definition B.2. For any ^)^ jajaj (sj) (j = 1, ■ ■ ■ , n) the operator 
TMVWicuOi) • • ■ 4>xL n <r n a n ( s n)) is defined by 

TMVWid! («l) ' ' ' ^x n a n a„ (s„)) 

:= Sgn(7r)^» x „(i)<7„(i)0».(i) ( S 7r(l)) 1 1 1 '0x 7r( „ ) <T^ ( „ ) a x( „ ) (*7r(n))) 

where 7r e 5„ is uniquely determined by the conditions that 

[-0x x(1) <r x(1) a, (1) (s T (i))] >T [V'x, (2) <T x(2) a x(2) (s 7r (2))] >T ••• >T [^x„ (n , <7„ (n) o„ (n) («7r (n) )] , 

and if there exist l\, l 2 e {1, • • • , n} with l\ < l 2 such that 

[ip^at^M] = bP* l2 a l2 a l2 (s l2 )},^ [^ jtTjaj {Sj)] (Vj € {/i + 1, l 2 ~ 1}) 

and 7r(m) = l\ with m G {1, • • • , n}, then 7r(m + 1) = l 2 . 

Using the ordering operator T\ we have the following expansion. 
Lemma B.3. For any t\,t 2 <E R with t\ < t 2 , 

e -(t 2 -ti)(H +Vx) 

= e -M)ff. +e -W.f (± f dsi...ds„Ti(T/ A (si)...T/ A (s„)) e * 1 ^. 
^i n! 

Remark B.4. Though the operator T\(V\(si) ■ ■ ■ V\(s n )) is defined only for 
(si,-- - ,s„) with Sj 7^ Sfe (j 7^ fc), we can consider Ti(Vx(si) • • • V\(s n )) as a 
Bochncr integrable function over [ii,i 2 ] n since the Lebesgue measure of the set 
{(si, ■ •• ,s„) e [*i,* 2 ]™ I 3j, k e {!,••• ,n} s.t. j 7^ fc and Sj = s fe } is zero. 
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Proof of Lemma \B.3[ Since the operator-valued function £ i— ► e (* 2 t i)( H o+4 y ^) j s 
analytic, we have 

It is sufficient to show that for all n G N and all £ G 

= (_l)» e -*2(Ho+en) f d Sl ...d S „T 1 (l/ A , 5 ( Sl )-.-l/ A , c ( S „))e* 1 ( ff « + ^), 

•/[ti,i2] n 

(B.2) 

where V\^(s) := e a< - Ho+ ^Vxe-< Ho+ ^ . In fact, substituting (|R2)l for £ = into 
(IB.ip gives the result. We show (|B.2I) by induction on n. 
By Lemma 12.31 we have 



e -(*2-ti)(iJo+ev x ) 



A e -(t2-*i)(Ho+fVx) = / dse -(l-s)(t 2 -ti)(H +^VA)(- ti _ t2 \pr Ae - S (t 2 -t 1 )(i?o+«V x ) 







ti 



which is (|B.2[) for n = 1. 

Let us assume that (|B.2[) is true for n — 1 (n > 2). 



e -(t2-ii)(/fo+Sn) 



= (-l)"- l (n- l)| e -*a(*o+£n) / dsids2 . . . dSn l 

JltiM"- 1 

■ l. 1 >*>...>.„_ 1 V A , e (a 1 )Vx,f(*2) • ' ' ^(s„_i)e tl ( ff «+« v ^ 
= (-l)"- 1 (n-l)! /""dsi P ds 2 --- /"" 2 ds„_ie- (t2 - Sl)(ffo+ ^V A 

. e -(si-S2)(i?o+«Vx)^ . . . e -(«„- 2 -«n-i)(H +eVx)y Ae -(s„-i-ti)(Ho+«V A )_ 



By writing £2 = soj ^i — s n and using Lemma 12.31 again we observe that 
'£ 



/so rs„-2 n 1 

d Sl ■ ■ ■ / d S „_! £ 

. e -( So - Sl )(H +£VA)y x . . . A L-isj-Si+^iHo+ZVx)^ y x 
. . . e -(s„-2-s„-i)(Ho+eVx)y A g-(s«-i-s„)(Ho+eVx) 

4G 



/•SO rSn-2 " _1 

(-l)»(n-l)!/ &!•••/ d S „_!^ 



( S o-si)(ffo+«V A )y A ... F' dse -( Sj - s )(H +tVx)y xe -(,s-s i+1 )(H +tVx)y x 

,-(s n -2-s n -i)(H +ZV x ) v ^ e -(s n - 1 - Sn )(Ho+tV x ) 
n-1 



" /'SO fSti-2 /-Sj 

(-l)"(n-l)!^ / dsi-../ ds„_i / dsl s 



rs re -2 

L s >si>--->s 

J - =0 "'s„ ^S n "'Sj + l 

. e -( So - Sl )(H +€Vx)y A . . . e -(sj-s)(ffo+«Vk)y Ae -(s-s; +1 )(ff +£Vx)^ 

. . . e -(sn-2-Sn-l)(Ho+eV x ) VAe -(s„_ 1 -S„)(Ho+eV x ) > 

Then by changing the index of the variables {sj, s \ j = 0, • • -n} we obtain 
= (-l)"7i!e- t2 ^ 0+ «^) / d S i---d S „l Sl >...> s „^(si)---^,c( S „)e tl(ffo+ «^ ) 

>/[ti,*a] n 

= (_l)« e -*>(Ho+CVO /" ds 1 ---ds n T 1 (V XA (s 1 )---V x , i (s n ))e t ^ Ha+ ^\ 

which completes the proof. □ 

Next we prepare some properties of the operators ipxo-( s ) an( i Vw(s)- Using the 
matrix {F(x<r, yr)}( x ( y r )grx{t j,} defined in (|2.2j) . we define the matrices F(a) 
(a = 0,1) by 

F if a = 0, 
-F* if o = l. 



F(a) :- 



Lemma [B.5I and Lemma [B.6I below follow [13 Lemma 3.2.1, Lemma 3.2.2]. How- 
ever, we give the proof to make this section self-contained. 

Lemma B.5. The following equalities hold. 

(i) For any (x, <r,a,s) ET x {f, j} x {0, 1} x M 

1*W(») = X; E e- sF W( XC r,yT)V y xa. 
yerre{T,i} 

(nj For any (x,cr,s), (y,r,<) e T x {|,J.} x K 

{iM^^M*)} - {^(*),^(«)} = 0, HUs)^yr(t)} = e^ F (yr^a). 
Proof. We see that for a € {0, 1} 

-^xaa(s) = e sHo (H ^aa ~ ^aaH )e~ sHo . (B.3) 

as 
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By using (IA.2j) we can show that for a £ {0, 1} 

-HqVW = - ^ ^ F(a)(x<7,yT)V>yTa +^xaa-ffo- (B.4) 

yerre{T,i} 

By combining (|B.3[) with (|B.4|) . we obtain a differential equation 

^~Vw(s) = - ^ X! ^( a )( xcr >y T )^yTa(s), 
S yer T e{T,|} 

for a G {0, 1}, which gives 

By using and (IA.2j) . the first equalities of fnj can be proved. Moreover, we 
see that 

{V'x ff («)»VVr(*)}= XI X e sF (xicri,XCT)e~* F (yr,X2CT 2 ){d ( T 1 ^ X 2 CT 2} 
xi,x 2 er C r 1 ,(72e{T,|} 

= E E e- tF (yr,x 1 a 1 ) e sF (x 1 a 1 ,x ( 7)-e( s -*) F ( y r,xa). 
xier<ne{T,4} 

□ 

For any linear operator A : F f {L 2 (T x {|,4};C)) -> F f {L 2 (T x {T,|};C)), let 
(A) Q denote Tr^-^^/Tre - ^ . For a set of the operators {ip xmaj (sj)}j=i> 

let Vxio-iai (si) • • • ip Xj ^ jaj (sj) ■ • • ip Xn a n a n (s n ) denote the product obtained by elim- 
inating V'xjo-.a, (sj) from the product ip Xl<Tiai (si) ■ • ■ ip Xn a n a n ( s »)- 

Lemma B.6. If n £ N is odd, /or any (xj , aj , aj , Sj) £ T x {f, i} x {0,1} x R 
(j = 1, • • • ,n) 

(Vitriol • • • V'x„ ( T„a„(s„)) = 0. 

If n £N is even, for any (xj, Oj, Oj, Sj) e T x {f, J,} x {0, 1} x K (j = 1, ■ • • , n) 

xi criai 

(Si) • • • ^x n <7 n a n (s n ))o 

n 

= E( _1 ) 3 ^ x i CT i a i( Sl )^^ a ^ S i))o(V'x2 CT 2a2(s2) ' ■•^x iff jOj(Sj) ' ■•^x„ CT „a„(s„)) - 

(B.5) 

Moreover, 

(V'xiffioi (51)^X2^202 ( s 2)> 

= E E (^ + e-' 3F(ai) )~ 1 (x 1 0- 1 , y T){Va 1 (si),V'x2 CT 2a2(s 2 )}. (R6) 

yerrefT,!} 
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e 



- 0Ho ^ ia M ■ ■■ip* n a n aM<t> G ML 2 ((T x {T, 1})';C)) 



(B.7) 



Proof. By using the ortho-normal basis Ur>f=o ^ m defined in (|A.1|) , we can write 
Tr^^ Vx lCTlQl V-x„ ff „a„ («„)) 

2L d 

= X ^ (0,e" /3ffo i/' Xl(Tiai (si)---V' Xii(Tiiaii (s„)0) F/ . 
Since for all s G R and m G {0, 1, ■ • • , 2L d } 

e sH ° (A m (L 2 ((r x {T,l}) m ;C))) c A m (L 2 ((r x {t,|}) m ;C)), 

we see that if n is odd, for any m G {0, 1, • • • , 2I/ d } and </> G B m 
or 



with / 7^ to, which implies that 

(0,e-' 3ffo Vx lCTl a 1 (si)---V'x„<T„a„(Sn)0) i , / =0. (B.8) 

The first statement follows from (IB.7I) and (IbT 



Let us assume that n is even. We see that 

X2CT2G2 
xi<ti ai X2CT2a2 

- ^x 2 <T 2 a 2 (s2)-0 xi (Ti ai (S3) • • • V'x„<r„a„(Sr ! ) 

xitTi ai X2CT2a2 (52)}^ X3CT 3 a3 n a n (^n) 

- V'x 2 <T 2 a 2 (s2){V' x 3 0-303 

+ V'x 2CT2 a 2 (s2)V' xi cti ai (S4) ■ • • V ; x„ C r„a„(s„) 

n 

= ^{V'xia-iQilsO^x.^a^SjOK-lF^aaa,^) • ■ • t/>x, aj aj (Sj ) • ■ • ^x„ CT „a» On) 
J=2 

(S2) ' ' "0x„ o -„a„(Sn)V l xi<ti ai 

which yields 

X2CT2a2 n a„ ( s ra))o 

+ (^x 2 <r 2 a 2 (s2) ' ' • V>x nO - ra a„(Sn)V , xi O - 1 ai(si)) 



= ^{^xttnoxOO.^^o^Si)} 

' (-l)- 7 ('0x 2 <T 2 a 2 (s2) ' ■ •V'x^j-ajOi) • • "0x„ O -„a„(Sn)) O - 



(B.9) 
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On the other hand, by Lemma [B.5l (|T|). 

^x™(s + /?) = E E e-™(xa,yr)^ yTQ ( S ). (B.10) 
yer T e{T,|} 

By using (jB.lOjl and the equality that Tr(AB) = Tr(BA) for any operators A, B, 
we observe that 

(V'x 2C r 2 a 2 (s2) ' ' "0x„a„a„(s„)V , x lO ' 1 a 1 (si)> o 

xi (Ji ai X.20-2CL2 (sa) ' ' ' V>x„ CT „a„(Sn))o 

= E E e_/3F(ai) ( X l cr l'y r )(V'yra 1 (si)V'x 2( T 2 Q 2 (s2) • • •V'x„ ct „q„(s„)) . 

yerre{T,i} 

(B.ll) 

By substituting (|B~TTj) into (|B~9| we obtain 

E E (kior-W + e-^Wx.yr)) 
yerre{T,l} 

(^2) ' ' 'V J x 7l (r n a ri (Sn))o 

n 

= E^ x i^ a i( s i)'^ x 3 ^ Q j( ;5 3 )J'(~ 1 )' : '^ x ^2a 2 (s2) ' ■■^ j a j a j {Sj) ■ ■ ■ <i/>x» CT „a» {s n )) - 
i=2 

(B.12) 

Let us define a unitary matrix M = (M(kr, xcr ))(k,T)er*x{T,J.},(x,CT)erx{T,J.} by 

M(kr,xa) :=^e-^). 

Then we have for all (k, r), (k, f ) € T* x {f, |} 

MFM*(kr, kf) = ~MF t M t (kT, kf) = 5 ki fe^r,f (B.13) 
where i?k is defined in (12. 8|) . The equality (|B.f 3j> implies that 
det(J + e"^) = det(J + e "/^™*) ^ 0, det(J + e^') = det(J + e ^' Mt ) ^ 0. 
Thus, for a = 0, 1 the matrix / + e~^ F ^ a > is invertible. The equality (|B.12|) leads to 

x 2 er 2 a 2 ( s 2) ' ' ' ^x„av,a„ ( s ri))o 

n 

= EE E (^ + '=^ F(ai) ) (xicn,yr){^ Tai ( Sl )^x iCT , a ,( Sj -)}(-lF 
i=2yerre{T,l} 

(S 2 ) • • • V^aj Oi) ' ' ' ^x„ CT „a„ (S„)) . 

(B.14) 

The equality (|B.6j) is (|B.14|I for n = 2. Then, by substituting (|B.6|) into (|B.14jl we 
obtain (TB31). □ 
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From now we show some lemmas involving the ordering operator T 2 . To simplify 
notations, let ipj denote ip Xjaja i (sj) for fixed variables (xj, <jj, a,j, Sj) € T x {|, J,} x 
{0,1} xR(j = !.•••»")• 
Lemma B.7. For any tt (E fSrtj 

r 2 (^i^2 ■■■i>n) = sgn(7r)T 2 (^(1)^(2) • ■ ■ VV(n))- (B.15) 

Proof. It is sufficient to show (|B.15|) for any transposition 7r as any permutation is 
a product of transpositions. Let us assume that 7r = (j,k), 1 < j < k < n. Let 
t, 77 € S„ be the unique permutations associated with the definitions of T 2 (^i ■ ■ - ipn) 
and T 2 (VV(i) ' ' ' respectively. 

T 2 (ipi ■ ■■ipn) = sgn(r)^ r (i) • • •VV(n). (B.16) 

12(^(1) ' ' '^(n)) = Sgn(?7)V' 7r (^(l)) • ••V'7r(»j(n))- ( B - 17 ) 

First consider the case that [ipj] 7^ [?/>&]■ Let A, B C {j + 1, • • • , k — 1} satisfy 
that [ipj] = [ipa] f° r an y a £ A, [ipk] = [^7] f° r an y 1 & B and [-0^] 7^ [t/> p ] for 
any p € {j + 1, • • • , k - 1}\A U B. 

It A,B 7^ 0, we can write A = {ai, • • • , a;}, 5 = {71, • • • , j m } with j + 1 < 
ai < • • • < ai < k — 1, j + 1 < 71 < • • • < 7 m < k — 1. By the definition of T 2 the 
product ipwfriti)) ' ' ' VVfrK")) is obtained by replacing ipjip ai • ■ • ip ai and i/> 7l • • • V> 7m V>fe 
by V'ai ' ' ' ^ai^j an d rfkipfi ' ' ' "07m respectively in the product VV(i) ' 1 ' ^r(n)- Thus, 
if we define cycles (1 , C2 G S„ by 

> _ ( j ai ■■■ a ; „i ai \ . _ / 7x 7 2 • • ■ 7™, fc \ 
1 V oix 0.2 ■■■ m j J' 2 \ k 71 ••■ 7 m _i 7 m /' 

the permutation r\ is written as 

V = K-Hibr. (B.18) 
On the other hand, Lemma IB. 51 (ii) ensures that 

"0a! ' ■■Ipailpj = (-l)VjV'ai '-'^an ^k1p 7l ' ■ ' ^ 7m = (-l)"V 7l ■ • ■ 1p lm 1p k . (B.19) 

By (pT6]l - (|BT9|l we see that 

T 2 (lp 7r (l) • • • V'Tr(n)) = Sgn(7T- 1 ClC2T)V' Cl ( C2 ( T (i))) • ' "0Ci(C2(r(n))) 

= (-1)!+'+™ sgn(r)(-l) i+m V- T(1) • • • ^ T(n) (B.20) 

= -T 2 (V>1 • • ■ V>n)- 

If A = or B = 0, by setting (i = Id and I = or ( 2 = Id and m = 0, 
respectively, we see that the equalities (|B.18|) and (|B.20[) hold true. 

Next consider the case that [ipj] = [4>k\- Let Ac {j + 1, • • • , fc — 1} be such that 
[ipj] = [tp q ] for any q e A and [^-] 7^ [t/^] for any 5 € {j + 1, • • • , k — 1}\A 

If A 7^ 0, we write A as A — {qi, ■ • • , q r } with j + 1 < qi < ■ ■ ■ < q r < k — 1. 
By the definition of T 2 the product VV(>j(i)) ' ' ' VVfjjfn)) is obtained by replacing 
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ipjip qi ■ ■ ■ ipq r ipk by ipkipqi ■ ■ ■ ipq r ipj in the product V>t(i) ■ ■ • ip T (n)- Thus, the permu- 
tation rj satisfies the equality 

77 = t. (B.21) 
By Lemma [B.5I (ii) the following equality holds. 

Ipklpqi ■ ■ ■ 1pq r 1p 3 = -Tpjlpqi ■ ■ ■ i>q T ^k- (B.22) 

By combining (|BTl6]l - (|BTf| with (|b721]) - ([BT22|) we have 

72(^(1) ' ' ' 1pn(n)) = Sgn(r)'0 7r ( T (i)) • • • -07r(r(n)) 

= - sgn(r)VV(i) • • • i>r(n) (B.23) 

= -T 2 (lpl ■ ■ -1p n ). 

By repeating the same argument as above without the term ip qi ■ ■ ■ %j) qr we can 
prove the equalities (|B .23|) for the case that A = 0, which completes the proof. □ 

Lemma B.8. Assume that n £ N is even and [ipi] >; [ipj] (Vj £ {2, 3, • • • , n}). The 
following equality holds. 

n 

(T 2 (V>! • • • ^„)) = ^(-l) j (T 2 (^^))o( T 2(^2 -Tr- i>n)) - (B.24) 

3=2 

Proof. For n = 2 the equality (|B.24j) is trivial. Assume that n > 4. Let t £ S n be 
the unique permutation associated with the definition of T^C^i ' ' • VVi)- 

T 2 (^l • --Ipn) = Sgn(r)VV(l) ■ • ' VV(n)- 

By assumption r(l) = 1. Moreover, Lemma IB .61 ensures that 

n 

(T 2 (V>1 • • • ^n)) = Sgn(r) ^(-l) 3 (V>lV>-r(i))o(VV(2) ■ ■ ■ VV(j) ■ ■ ■ VV(n)) - (B.25) 

Let us fix j £ {2, 3, • • • , n}. Let it £ S n be such that 

Ml), 7r(2),7r(3), • • • , 7r(n)) = (1, r(j), r(2), • • • , • • • , r(n)), (B.26) 

where 't(j')' stands for the omission of the number t(j) from the row 
(t(2),t(3),- • • , r(n)). Then, we have 

sgn(7r) = (-l)^- 2 sgn(r) = (-1)^ sgn(r). (B.27) 

On the other hand, we can write {1, ■ • • , t(J)} = {h, ■ ■ ■ ,^-2} with 

2 < h < h < ■ ■ ■ < ln-2 < n. There exists r\ £ 5„_ 2 such that 

(l v (i),l v (2),- ■ ■ J v (n-2)) = (t(2),t(3),--- ,t(j),-- - >r(n)). (B.28) 
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(B.31) 



By (|R26) and (fB~28l) we obtain 

(7r(l),7r(2),7r(3), ■ • ■ ,7r(n)) = (1, r(j), l ri{1) , ^ (2) , ■ ■ ■ ,l v ( n -2)), 
which implies that 

sgn(vr) = (-iy(i)- 2 sg n(r ? ) = (-l) r « sgn(??). (B.29) 
By (|B~27l) and (|B~29f we have 

sgn(r) = (-1) T « sgn(ry). (B.30) 

Note the equalities that 

(T 2 (V>iVv(j)))o = (V'lV , rO')) 0! 

(T 2 (^ 2 ' ' ' VV(j) ' ' ' tpn)) = Sgn(7 ? )(V' T (2) ' ' ' 1p T (j) ■ ■ ■ VV(n)) - 

By substituting (|B.30|I and (|B.31| into (|B.25jl we see that 

n 

(T 2 (Vi • • • V>n)) = ^(-l) T W(T 2 (^ 1 ^ a) ))o(T 2 (^ • • • Vv^ • • • lM) 

n 

= £(-l) J '<T a (V>iVi))o<Ta(^ • • • • • • V-n)> , 

which is JH22). □ 
Lemma B.9. For a/Z Xj, S T, (7j,Tj G {t, 4}. Sj, tj el (j = 1,2, • • ■ , n), 

( T 2(V£ 1CTl (si>i/Vi-ri(*l) ' ' ' V4„ ( r„( s n)V'y„r„(in))>0 , R 

= det((T 2 « J ^(« i )^(<fc))>o)i<M<«- ' 

Proof. We show (|B.32I) by induction on n. The equality (|B.32[) is obviously true 
when n = 1. Let us assume that (|B~32|t is true for n - 1 (n > 2). 
Lemma IB. 71 implies that for all it £ S n 



(?2W'x 1 <7i( a l)V'ym (*0- ••V'x„<T„(Sn)V'y„T„(in)))o 

= ( T 2(C (1)( r, (1) ( s 7 r(l))'0y^ (1) T, (1) (^(1)) ' ' ' ^x w(n) ^(„) ( S 7r(n))V'y T („)T T („ ) (*7r(n)))) 

= (-i)"(T 2 (^ cl) ^ (1) (^ (1) )C (I) ^ (I) (Mi)) 

• • • i>y*( n )T„ M (*ir(n))^„ (n) tr, (n , ( S 7r(n)))) ; 

(B.33) 

and 

dct ( (T 2 (V£ (Sj)lpy kTh (tk))) )l<j,k<n 

= det((T 2 (^ (j)(T7r0) (s 7r ( 3 ))V'y, (fc) r 7r(fc) (^(fc)))>o)i<i,fc<n (B.34) 
= (-1)" det((T 2 (^ y7r(fc)T7r(fc) (^(fc))C (3) <r x(3) ( s ^(i)))>o)i<i^<«- 
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The equalities (|B.33|) and ()B.34[) enable us to assume that 

h^ lCTl (si)] b bPZjajlsj)}, [i> yjT /(tj)} (Vj G {1,- • • ,n}) without losing generality in 
the following argument. 

By using Lemma IB. 71 Lemma IB.81 the hypothesis of induction, and the fact 
that (V>xaWVCv(i')>o = °> we have 

(^(d^CsOV'ym (h)- ' •C l( r„( s n)V'y„T„(i™))) 

n 
3=1 

• (7 1 2(V'y 1 r 1 (il)C 2 <T 2 ( s 2)V'y2r 2 (i2) ' ' ' V4 3 <x 3 (SjOVVjT; Oj) ' ' ' ^x„<r„ ( s n)VVr,r„ (*n))) 
n 

= E(- 1 ) i ~ 1 < T2 ^>0VW^ 

= det((T 2 (V'K 3 . (T .(s i )^y fc T fc (tfc)))o)i<i,fe<n, 
which concludes the proof. □ 
Lemma B.10. For all x, y G T, a, r G {t, I}, ijeM 

( T 2('0x<T(a ; )V'yT(y))> o = C{xo-x,yTy), (B.35) 
where C(xcrx,yrj/) is defined in (|2.7|) . 
Proof. By the dehnition of T2 , we have 

(^(VCO^yr^iOo = (V'x ( T(2 ; )^yr(2/)) lx-y>0- (^(^Co-O^O^-jKO- ( B -36) 

By Lemma lB31(|n|l. (fB~6) and (jB~T3]) . we see that 
(V4t<>)VVt(2/))o 

= E E (/ + e^ t ) _1 (x ( 7,xV'){Cv'W»V'yr(y)} 

x'erve{t4} 

= E E (^ + e /5Ft ) _1 (xa,xV) e ^ F (yr,xV) 
x'erve{t,|} 



L d °k,k L j + e /3.E £ Ld 2^ 1 + e /3£ k ' 

k,ker* k 6 r * 

(B.37) 
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= E E (/ + e-H~V,yV){^ y v<(y),C>)} 

y'er T 'G{T,l} 

= E E (/ + e-^)" 1 (yr,yV) e M^(yV,x ( 7) 
y'er T 'e{T,l} 

= ({l + e-P F y X e^ F \ (yr,xa) 

2^ °k,k e e , , .-/JJS. " r d 



k,ker* k er* 

(B.38) 

By combining (|R37|) and (|BT38|) with (|B~36)) . we obtain (|B~35f . □ 

We have prepared all the lemmas necessary to prove Proposition 12.41 
Proof of Proposition\2.4\ By applying Lemma |B. 31 for ti = 0, t% = j3 we have 



+ e -wf f ds 1 ---ds n T 1 (V x (s 1 )---V x (s n )) 

CO . 

e + e -P B °Y j (-l) n / dsi---ds„l sl >...> Sn FA(si)---F A (s„). 

n=l •'fO-fll" 



[0,/J]" 

(B.39) 



By (|B.39|) , Lemma IB. 51 (fn) , the definition of T2 and Lemma IB. 71 we see that 

Tr e- 0Hx 
Tr e-P^ 



" 1 + §n - E 



• (V'; iT ( s i)V'y 1 i('Si)V'w 1 i(si)i/'z 1 T(si) ' ' •C„T( s ») ? /'y re l( s ™) 1 /'w„i(s„)V'z„T(sn))o 

= 1 + E 3 3 i r E f dX 2j -l U X2] _ u k 23 ,y 2j _ ! ,y 2j 

lxi>a3>"->X2n-i 

(-l) n (V>; iT (a;i)?/;; 2i (a;i)-!/; yiT (a;i)7/; y2i (a;i) ■ • • 

' C 2 „_ 1 T( a; 2™-l)C 2 „x( a; 2r l -l)^y 2 „-iT( a; 2 n -l)^y 2 „l(a;2«-l)>o 
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n=l J y x 2 j-i,x 2 j,y2j-i,y2jer 

•V';2„-lT( a;2 «- 1 )^y2™-iT( a; 2n-l)C 2 „l(2;2n-l)V'y2„l(2 ; 2n-l)))o 
1 n 



i + > -: n 



n=l 



/3 



^ ] ^ ] / ^ 2 '2j-l^cr2j_l.T^O'2j,l^X2j-l,X2j,y23-l.,y2j 

X2j-i,X2j,y23-i,y2 3 -er a 2 j-i,<T2j£{1d} 
(^(V'xitrx (*l)VWi (^l)C 2 a 2 ( x 2)fpy 2 
■ C 2 „_ lCT2 „_ 1 (^n-O^ysn-KTsn-i (a;2„-l)C 2n cr 2n ( a; 2n)'0y2„<T2„ («2n)))o 



a: 23= :c 2j-l 
Vj€{l, -- 



Then by using Lemma IB. 91 and Lemma IB. 101 we obtain the series (|2.6j) . □ 

C Diagonalization of the covariance matrix 

In this part of Appendices we diagonalize the covariance matrix 
(C/^xcrir, y-ri/))( x ax ^ .(y,T,i/)erx{T,l}x[o./3)h an< ^ calculate its determinant. The fact 
that the determinant of the covariance matrix is non-zero, which is to be proved 
in Proposition IC.71 verifies the well-posedness of the Grassmann Gaussian integral 
defined in Definition 13.51 

For convenience of calculation we assume that h G 2N//3. Define the sets Wh 
and Mh by 



nil u.- -nil [ . Xf h := <J W€ ^(2Z+1) 



irh < uj < Tch 



Note that \W h = 2j3h and \M h = (3h. The assumption that h G 2N//3 ensures the 
equality 

M h = W h \^. (C.l) 

The set Mh is seen as a set of the Matsubara frequencies with cut-off. 
For / G i 2 ([-/3, p) h ; C) we define / G L 2 (W^.; C) by 

/H : =^ E e-^VW- 
te[-p,P) h 

Lemma C.l. For any f G L 2 ([-/3, C) 

• f W = ^ E Vt €[-/?,/?)*. 



5G 



Proof. If t = -f3 + s/h with s G {0, ■ ■ ■ , 2/3/i - 1}, 

i £ e^/M = ^ E E e^e-^fiu) 

' u<£W h f uj£W h ue[-p.j3) h 

2/S/i-l 20/1-1 



- —-.'II J. J. / , 

1 ^7 ^ e i(-Kh+-Km/P){a/h-l/h) f ( _p + i 



2Bh ^ ^ J \ ' h 

2/3fc-l 2/3/1-1 . x 

— E e _i7r ^ - ^ E e 47rm(s " i)/(/3,l) / ( -/?+ - ) 

z=o m=0 ^ ' 



2/3/t-l 



Z=0 



E e-***-^.,,/ ( "0 + ~ J = / {-P + J) = /(*) 



□ 



Lemma C.2. Iff G L 2 ([-[3, 0) h ; C) satisfies /(*) = -f{t + [3) for all t G [-/3,0)h 
with t < 0, 

= ^ E e ^/»> vt e h*./ 3 )*- ( c - 2 ) 

Proof. Take any o> G fl 27rZ//3. By assumption we see that 

K") = { E e_ia "/(i) + ^ E e ~^/w 

te[-0,0) h \[o,/3) h te[o,p) h 

= -\ E e " Mt + \ E e ^ 4 /w (c.3) 

te[-/3,/3)h\[0,/3) h te[o,/3) fc 

= E e"*-(*-«/(*) + ^ E e-*-*/(*) = 0. 
te[o,/3) h te[o,/3) h 

Then, by l|CTT]) . (|C~3|) and Lemma EH] we obtain ljCT2|l . □ 

Let us define g k G L 2 ([-/3, C) (keP) by 

tE k J 1*>0 1*<0 



1 + e /3Sk 1 + e -PE k 



Note that the function g^ satisfies the anti-periodic property gk(t) = —gk(t + (3) for 
all t G [-0,0)h with i < 0. 

Lemma C.3. For all t G [— /3,/3)h 

9±(t) = p E wi-rWW)' (C - 4) 
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Proof. By Lemma [C. 2 1 



= Ta E (c.5) 



Moreover, we observe that for ui £ Mh 



i . . P tE k -i „tEk 

- /, ,\ _ 1 --tot e I „-wt e 

9kM--- ^ e _ pEk+ - ^ e TT ^- £k 

t6[-/3>/3)h\[0,/3) h *e[0,/8) h 



-iuit ° 



_A V P -^(«-^) 6 i V 

te[o,/3)h te[o,/3) h 

£ \ p-«w*__r__ = \ t(-iu+E k ) 

h ^ l + e P E ± ftfl + e^*) Z- 



e 



2 



(C.6) 



fl(\ — e -iu/h+E k /h\ ' 

The equality ([04)) follows from |C3)) and (fUT6|) . □ 



By substituting the characterization of (?k given in Lemma IC.3I into 

C h ( X ax,yry) = ^ £ e l < k ^ 5k (x - y) 
ker* 

we obtain 

Lemma C.4. For any (x, a, x), (y, r, y) e T x {|, J,} x [0, /3)fe, 

C^x^yry) = ^ ^ ^ ft(l - e-W^AV ( CJ ) 
In order to diagonalize C/j, we define a matrix 

Ytkrcu^ax) := feg P »(k,x) p -»^ 
Lemma C.5. TTie matrix Y is unitary. 

Proof. Assume that w = — %h + it/ (3 + 2wm/f3, Cj — —irh + ir/[3 + 2irrh/ (3 with 
m, rh £ {0, 1, • • • , 8h — 1}. Then we observe that 

xerze[0,/3) h 



(3h f-^ 
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Let x = s/h, x — s/h with s, s € {0, 1, • • • , j3h — 1}. 



Y*Y(xax,xax) = J2 

ker* ugM, 



i(k,x— x) itu(x-x) 

(3hL d 



e 

0h-l 



^<T,g-^x,x i{—Kh+-K/0+2-Km/P){s/h-s/h) 
H m=0 

0h-l 

Og.o-Ox.x r i(-irh+Tr/0)(a/h-S/h) e i2Tvm(s~s) i »'/< i 



m=0 

A -A -A . P i (- 7r ' i + Tr /0)( s /' l -s/' 1 ) — A . A -A - 



By using the matrix Y and (|C7[) we can diagonalize Ch as follows. 
Lemma C.6. For all (k, r, w), (k, f , u) € T* x {f, 4} x M h) 

(YC h Y*) (k™,kf£) = 6 T , f 6 k £6 u ,z i - e _ J /h+gk/h ■ 

Finally we calculate the determinant of the covariance matrix C/,.. 
Proposition C.7. For am/ /i <E 2N//3 

n ker .(i + e^) 2 

Proof. Since { e ~^/h+E k /h | ^ g M/ j ig the get of &n the ^ hth roots of _ e /3is kj 

z 0h + ( ,0E k = n ^ _ 

LueM h 



□ 



for all z G C. Especially, 



n (1 - e ~ lul/h+Ek/h ) = 1 + e f3Ek . (C.8) 



By Lemma TC. 5 1 Lemma IC . 6 1 and (|C.8[) . we see that 
detC h = det(YC h Y*)= n n n 



ker* ae{ T ,l} «eM k 1 - e -^/ fe + E */>> riker* (1 + e^*) 2 ' 

□ 
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